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ABSTRACT 

A  new  approach  to  defining  inlet  distortion  factors  for  axial -flow 
compressors  is  presented.  A  summary  of  past  analyses  of  inlet  distortion, 
along  with  a  brief  review  of  flow-distortion  definition  and  compressor 
performance,  provides  a  background  for  the  proposed  theory.  Inlet  flow- 
distortion  in  terms  of  total -pressure  fluctuations  at  the  compressor  face 
is  converted  to  vorticity.  The  effects  of  vorticity  on  the  internal 
aerodynamics  of  the  compressor  are  investigated  with  some  approximate 
calculations  included  to  indicate  the  validity  of  this  approach.  Refine- 
ments and  alternate  solutions  to  the  theory  which  could  lead  to  definition 
of  a  more  comprehensive  and  reliable  definition  of  stall-inducing  inlet 
distortion  are  included. 
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I.   INTRODUCTION 

Engine  stall  has  always  been  a  problem  of  varying  magnitude  in  jet 
powered  aircraft  employing  axial-flow  compressors.  At  present  aircraft 
gas  turbine  performance  is  largely  limited  by  flow  separation  and  surge 
in  the  compressor.  Generally,  highest  compressor  efficiencies  occur 
near  the  point  of  flow  breakdown  (stall),  and  so  the  design  margin 
between  normal  operation  and  stall  is  kept  to  a  practical  minimum.  A 
compressor's  ability  to  remain  below  this  stall  limit  is  significantly 
affected  by  the  amount  and  degree  of  distortion  in  its  inlet  flow. 
Inlet  flow  distortion  has  been  called  by  many  names,  all  of  which  basi- 
cally describe  variations  in  total-pressure,  total -temperature,  or 
velocity  over  the  inlet  cross-section.  Previously,  inlet  flow  distor- 
tion has  been  at  a  level  which  the  compressor  could  handle  without  a 
serious  degradation  in  performance.  However,  with  the  greater  speeds 
and  performance  demanded  of  today's  high-specific-energy  turbomachines, 
inlet  distortion  at  the  compressor  face  can  alone  be  great  enough  to 
cause  compressor  stall  and  surge.  One  prominent  source  of  distortion 
has  been  traced  to  shock-wave/boundary-layer  interactions  within  the 
inlet.  Others  are  steam  ingestion,  pressure  transients  due  to  afterburner 
light-off,  ordnance  discharges,  etc.  The  magnitude  of  the  problem  in 
some  recent  aircraft  designs  has  caused  much  research  to  be  conducted  in 
this  area,  specifically  correlating  maps  of  total-pressure  variations  at 
the  compressor  face  to  engine  surge.  The  results  of  this  work  have 
yielded  an  immense  amount  of  empirical  data.  However,  reduction  of  these 
data  has  yet  to  provide  a  technique  which  adequately  describes  the  distor- 
tion accommodation  of  a  given  engine-inlet  configuration  over  its  entire 
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range  of  operation.  In  contrast  to  empirical  correlations  yielding  a 
statistical  yes-no  answer  to  the  question  of  whether  or  not  the  engine 
will  stall,  methods  are  needed  which  link  stall  criteria  for  a  blade 
element  to  distortion  patterns.  This  paper  presents  a  new  approach  to 
analyzing  already  existing  data  and  provides  evidence  of  a  better 
insight  into  the  actual  fluid  mechanics  involved  in  the  problem. 
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II.     STATEMENT  OF  PROBLEM 

A.     AXIAL  COMPRESSOR  PERFORMANCE1 

The  performance  of  an  axial -flow  compressor  such  as  found  on  air- 
craft turbojet  engines  is  best  described  by  its  steady-state  performance 
map.     Figure  1   illustrates  a  typical   compressor  map.     "Pressure  ratio" 
is  exit  total -pressure  to  inlet  total-pressure  and  "corrected  airflow" 
is  engine  mass-flow  rate  corrected  for  total -temperature  and  total- 
pressure. 
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FIG.  I        COMPRESSOR     PERFORMANCE     MAP 
J85-  GE-  13  [adapted  from  Ref.2J 


^The  information  and  examples  in  this  section  were  adapted  from  a 
report  by  Brimelow  [Ref.  1]. 
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The  surge  line  defines  the  maximum  pressure  ratios  which  the  compressor 
can  achieve  throughout  its  range  of  rotor  speeds.  Operation  above  this 
line  is  accompanied  by  a  drastic  reduction  in  pressure  ratio.  This  is 
called  "surge,"  and  it  is  caused  by  excessive  aerodynamic  blade  stall. 

The  stability  of  a  system  is  determined  by  the  slope  of  its  pressure 
ratio/flow  characteristic.  Consider  steady  operation  at  point  M-j  of 
Fig.  2. 
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COMPRESSOR  /  RECEIVER 
MATCHING       CHARACTERISTICS 


Now  if  the  flow  became  unsteady  and  a  perturbation  were  to  cause  a 
reduction  in  airflow  to  point  S-j  (constant  speed),  the  compressor  would 
be  delivering  an  instantaneous  pressure  greater  than  that  in  the 
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receiver.  Effectively  the  compressor  back-pressure  has  been  lowered, 
and  the  flow  would  increase  until  a  match  is  again  reached.  Similarly 
for  an  increase  in  flow  to  point  S2,  the  increased  back-pressure  would 
cause  a  return  to  stable  operation  at  point  M, .  However,  if  the  compres- 
sor is  matched  on  the  positive  slope  of  the  characteristic,  by  the  same 
reasoning,  unsteady  conditions  exist  and  perturbations  are  aggravated 
to  "surge." 

A  stall  margin,  defined  as  the  vertical  distance  between  the  nominal 
(design)  operating  line  and  the  surge  line,  can  be  thought  of  as  a 
measure  of  the  compressor's  ability  to  accommodate  flow  transients  which 
instantaneously  load  the  machine.  This  margin  is  dictated  by  overall 
engine-inlet  design  and  is  affected  by  many  different  phenomena.  Unfor- 
tunately the  point  of  maximum  efficiency  in  a  compressor  is  just  prior 
to  stall  or  just  below  the  surge  line  on  the  performance  map,  and  there- 
fore a  trade-off  between  efficiency  and  surge  margin  is  required.  The 
expected  flow  transients  which  the  compressor  must  be  able  to  accommodate 
will  affect  the  design  efficiency. 

Flow  transients  initiated  within  the  engine  system  are  normally 
caused  by  power  level  (throttle)  movement.  Transients  can  also  be 
initiated  by  flow  distortion  at  the  compressor  face.  Of  the  various 
flow  distortion  parameters,  Brimelow  [Ref.  l]  states  that  inlet  total- 
pressure  and  total-temperature  have  been  identified  as  having  the  greatest 
effect  on  engine  stability,  which  is  to  a  large  extent  the  compressor's 
ability  to  remain  below  the  surge  line.  Total -pressure  variations  can 
be  steady-state  spatial  distortion,  time-varying  spatial  distortion, 
in-phase  pulsations  which  are  spatially  uniform  over  the  compressor  face, 
or  a  combination  of  these.  Low  frequency  spatially  uniform  (planar) 


14 


pressure  fluctuations  (on  the  order  of  20  CPS)  do  not  affect  the  com- 
pressor but  do  have  a  detrimental  effect  on  the  inlet  and  engine  control 
systems.  High  frequency  planar  fluctuations  (more  than  one  per  revolu- 
tion) affect  compressor  stability  much  like  time-varying  spatial 
distortion. 

Spatial  distortions  can  be  further  segregated  into  radial  and  cir- 
cumferential inlet  flow  distortion.  This  separation  into  circumferential 
and  radial  distortion  was  arbitrary,  somewhat  motivated  by  geometric 
simplicity.  A  radial  total -pressure  distortion  such  as  shown  in  Fig. 
3(a)  produces  yery   little  variation  in  inlet  static  pressure.  The  net 
result  is  a  low  velocity  in  the  low  pressure  region  and  a  high  velocity 
in  the  high  pressure  region. 
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RADIAL     PRESSURE      DISTORTION 
[adapted  from  Ref.  l] 
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In  the  low  velocity  region  the  blades  will  operate  at  higher  incidences 
and  higher  total-pressure  ratios.  The  net  result  will  be  a  reduction  in 
the  compressor  stall  margin. 

For  the  circumferential  type  distortion  shown  in  Fig.  4(a),  the 
associated  surge  line  reduction  is  a  function  of  the  basic  stage-match 
of  the  compressor,  the  axial  spacing  between  blade  and  vane  rows,  and 
the  unsteady  response  of  the  rotor  blades.  Since  the  blade  rows  restrict 
circumferential  redistribution  of  flow,  each  circumferential  segment 
performs  essentially  independently. 
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All  segments  are  required  to  deliver  the  same  exit  static  pressure  so 
that  in  the  area  of  low  inlet  pressure  there  is  a  high  pressure  ratio 
and  vice-versa.  Correspondingly,  the  high  pressure  ratio  forces  a  low 
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axial  velocity  in  this  area  and  a  high  velocity  in  the  segment  of  high 
inlet  pressure.  These  variations  in  velocity,  produced  less  than  one 
duct-diameter  upstream  of  the  compressor,  result  in  static  pressure 
fluctuations  at  the  compressor  face.  This  mechanism  of  varying  veloci- 
ties (incidences)  can  induce  surge  when  the  average  pressure  ratio  of 
the  distorted  flow  is  less  than  the  surge  pressure  ratio  for  undistorted 
inlet  flow.  The  effects  of  circumferential  pressure  pulses  on  a  TF-30 
engine  are  presented  in  Ref.  3.  It  was  found  here  that  the  amplitude  of 
a  pulse  necessary  to  stall  the  engine  varied  inversely  with  pulse  dura- 
tion and  that  the  engine  was  most  sensitive  to  distorted  sectors  of  180 
to  240  degrees. 

The  distortions  which  are  encountered  in  actual  compressor  applica- 
tions such  as  jet  aircraft  inlets  are  neither  radial  nor  circumferential 
but  a  combination  of  the  two.  Furthermore  the  magnitude,  location, 
shape,  and  duration  of  these  pressure  variations  change  continually. 
It  is  virtually  impossible  to  evaluate  precisely  the  effect  that  these 
distortions  have  on  engine  stability  because  there  are  so  many  distor- 
tion patterns  significant  to  each  compressor  surge.  Therefore  it  has 
been  necessary  through  experimental  programs  to  develop  generalized 
correlations  between  pattern  characteristics  and  corresponding  effects 
on  compressor  stability.  Correlating  the  shape,  location,  spatial 
extent,  and  magnitude  of  the  inlet  distortion  requires  analyzing  many 
radial  and  many  circumferential  patterns  and  relating  the  overall  picture 
to  engine  stability.  From  this  is  generated  a  distortion  factor  which  is 
a  measure  of  the  associated  loss  in  stability-margin  for  that  given 
propulsion  system.  Because  compressor  geometry  and  blading  are  not 
considered  in  arriving  at  this  parameter,  the  result  cannot  be  applied 
to  other  systems. 
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For  a  more  detailed  and  complete  analysis  of  this  subject,  consult 
Brimelow's  excellent  work  in  Ref.  1. 

B.  PRESENT  STATE  OF  INLET  DISTORTION  ANALYSES 

Until  a  few  years  ago,  engine-inlet  compatibility  limits  were  success- 
fully determined  by  evaluating  an  engine's  tolerance  to  steady-state 
distortion.  Normally  this  was  done  with  screens  or  similar  obstructions 
placed  in  the  inlet  and  the  compressor  face  instrumented  with  low-response 
pressure  probes.  The  engine  was  operated  under  all  conditions  of  interest, 
and  by  comparing  pressure  data  at  the  compressor  face  with  the  engine's 
performance,  distortion  factors  were  derived.  It  was  known  that  unsteady 
flow  was  being  produced  by  these  test  devices,  but  no  attempt  was  made  to 
evaluate  its  effect.  However,  when  this  technique  was  applied  to  the 
TF-30  turbofan  engine  during  the  F-111A  flight  test  program,  the  data  did 
not  correlate  well  nor  was  the  method  found  valid  for  other  fan  type 
engines  at  that  time.  An  improved  steady-state  distortion  factor,  K .  ,    ' 
was  developed,  but  it  too  fell  short  of  providing  the  needed  total  evalu- 
ation of  inlet  distortion.  It  was  becoming  more  evident  that  there  was 
another  significant  variable  involved,  namely  turbulence.  Turbulence  can 
be  described  as  flow  with  random  pressure  fluctuations  and  can  be  quanti- 
fied as  a  fluctuation  in  RMS  total-pressure,  (aPdms^^pt^  Wlt'1 
increased  aircraft  speeds,  turbulence  has  taken  on  greater  importance 
since  there  is  a  significant  amount  of  turbulent  flow  associated  with 
the  boundary-layer/shock-wave  interaction  within  supersonic  inlets.  It 
was  felt  that  steady-state  distortion  indexes  were  still  of  value,  but 
a  correlation  with  turbulence  was  necessary.  One  such  empirical  relation 
showing  a  decrease  in  distortion  tolerance,  K,2,  with  increasing  turbulent 
energy  is  depicted  in  Fig.  5. 
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LOSS     IN     STEADY- STATE      DISTORTION 
TOLERANCE      WITH    TURBULENCE 
[adapted    from     Ref.    4 

Turbulence  has  two  primary  effects  on  axial  compressor  performance. 
First,  turbulent  inlet  flow  causes  a  reduction  in  compressor  surge  line, 
its  sensitivity  being  a  function  of  turbulent  frequency  and  amplitude. 
Second,  a  reduction  in  flow  capacity  accompanies  inlet  turbulence.  These 
two  phenomena  were  observed  in  tests  performed  by  Brimelow  and  Plourde 
[Ref.  5]  on  a  TF-30  turbofan  engine.  This  all  led  to  the  development  of 
new  distortion  indexes  such  as  the  Kqa  and  IC  parameters  of  Ref.  6. 
However,  lack  of  correlation  caused  Burcham  and  Hughes  to  raise  doubt  as 
to  the  usefulness  of  the  conventional  turbulence  definition,  (aPrmq)/(Pt)- 
Bellman  and  Hughes  [Ref.  7]  found  no  consistent  relationship  at  all 
between  turbulence  factors  and  surge,  and  in  some  cases  surge  was  experi- 
enced as  turbulence  level  was  decreasing.  Therefore  it  was  concluded  here 
that  peak-to-peak  pressure  fluctuations  at  key  areas  of  the  compressor 
face  are  more  important  in  surge  prediction  than  are  turbulence  levels. 
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Also  extensive  experimental  data  have  been  obtained  from  a  J85-GE-13 
turbojet  placed  in  a  supersonic  wind-tunnel  running  at  Mach  numbers 
close  to  2.5.  By  varying  the  characteristics  of  the  engine  inlet  duct, 
a  desired  level  of  time  variant  flow  distortion  was  produced  through 
internal  shock  interactions  [Ref.  2].  Calogeras,  Burstadt  and  Coltrin 
found  for  this  J85  configuration  that  time-averaged  distortion  could  be 
well  below  a  stall-inducing  value  while  instantaneous  peaks  exceeded 
this  value.  In  addition,  severe  pressure  fluctuations  could  result  in 
an  instantaneous  distortion  of  a  type  quite  different  from  the  time- 
averaged  distortion. 

Of  these  few  representative  programs  discussed,  there  are  two 
striking  similarities  worthy  of  note:  (1)  Each  is  dedicated  to  using 
inlet  flow  distortion  described  by  raw  total -pressure  data  to  predict 
compressor  surge.  No  attempt  at  determining  the  cause  or  nature  of  the 
stall  is  made.  (2)  There  is  much  uncertainty  and  disagreement  associated 
with  the  present  methods  of  predicting  distortion  induced  stall.  A 
universal  distortion  parameter  has  not  been  agreed  upon. 

Obviously  a  better  understanding  of  the  internal  fluid  dynamics, 
while  a  difficult  task  at  best,  would  nonetheless  allow  a  more  compre- 
hensive generation  of  a  stall-inducing  distortion  parameter.  Total- 
pressure  is  the  easiest  and  perhaps  the  only  flow  parameter  which  can  be 
measured  at  the  compressor  face  with  sufficient  accuracy  and  frequency 
response  to  gain  meaningful  data.  The  vast  amount  of  pressure  data 
available  due  to  the  previously  mentioned  programs  provides  an  exten- 
sive empirical  data  base  with  which  to  treat  this  problem.  In 
evaluating  these  data,  inclusion  of  time  and  temperature  dependence 
could  only  improve  the  usefulness  of  the  results.  Therefore  it  is 
proposed  that  vorticity  become  the  basis  of  a  "universal  distortion 
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factor."  This  thesis  will  show  how  existing  total-pressure  maps  along 
with  total -temperature  and  time  data  can  be  transformed  into  vorticity 
maps.  Analysis  of  the  resulting  vorticity  maps  will  provide  evidence 
as  to  the  flexibility  of  this  method. 

Vorticity  is  more  closely  related  to  fluid  dynamics  than  is  total 
pressure  alone.  Compressor  stall  and  engine  surge  are  flow-separation 
phenomena;  therefore  the  key  to  understanding  their  causes  lies  in  the 
internal  dynamical  situation. 

C.  VORTICITY  APPROACH 
1 .  Compressor  Stall 

The  most  significant  phenomenon  affecting  the  performance  of  an 
axial-flow  compressor  is  stall.  If  the  flow  rate  through  a  compressor 
blade  row  is  decreased  while  the  rotor  speed  is  maintained  constant, 
angles  of  attack  on  the  blades  will  increase  and  eventually  flow  sepa- 
ration will  occur.  This  condition  is  called  stall  and  is  similar  in 
some  respects  to  the  stall  of  an  isolated  airfoil;  differences  are  due 
to  the  adjacent  blading  in  a  cascade.  Stalled  operation  of  a  compressor 
is  accompanied  by  a  drastic  decrease  in  efficiency  and  excessive  oscil- 
lating blade  loads  capable  of  causing  structural  failure.  Stall  occurs 
in  the  range  of  unstable  compressor  operation  and  hence  can  trigger 
surge.  Usually  a  compressor  stalls  in  an  asymmetric  manner;  that  is, 
well  defined  regions  of  the  annul  us  are  stalled  while  the  remaining  area 
is  unstalled;  see  Fig.  6.  The  flow  through  these  stalled  patches  is 
severely  restricted  in  accordance  with  the  decrease  in  total  flow-rate 
through  the  compressor,  thus  the  flow  is  diverted  around  this  blockage. 
Effectively  this  increases  blade  angle-of-attack  on  one  side  of  the 
restriction  while  decreasing  it  on  the  other  side.  In  Fig.  7  blade  No.  1 
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FIG.  6 


TYPICAL      STALL       PATTERNS 


Decreasing 
Angle-  of  -  Attack 


Increasing 
Angle-  of  -Attack 


FIG.  7         STALL    CELL     GROWTH    AND 
ROTATION    IN    A    BLADE    ROW 
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will  soon  become  unstalled  while  at  blade  No.  5  flow  separation  will 
commence.  This  condition  is  termed  "rotating  stall"  since  it  propagates 
circumferentially  in  a  direction  opposite  to  rotor  rotation.  In  an 
absolute  frame  of  reference  this  speed  of  circumferential  propagation 
is  30-60  per  cent  of  the  rotor  speed  [Ref.  8].  When  these  conditions 
become  severe  enough  to  cause  net  flow-rate  fluctuations  with  time,  the 
situation  is  called  "surge."  This  describes  operation  in  the  unsteady 
area  of  the  compressor  performance  map  in  Fig.  1.  Valensi  [Ref.  9]  has 
made  a  thorough  investigation  into  these  phenomena  with  some  excellent 
photographic  results.  Of  the  various  analytical  methods  available  for 
predicting  compressor  stall,  the  most  generally  accepted  employs  the 
limiting-blade-loading  criteria  of  Lieblein,  Schwenk,  and  Broderick 
[Ref.  10].  Their  "limiting  diffusion  factor,"  D,  is  a  function  of  the 
velocity  relative  to  the  blade  and  the  circulation  about  the  element. 
If  a  constant  axial  velocity  assumption  is  made,  this  simplifies  to  a 
function  involving  the  flow  turning  angle  and  cascade  solidity.  For 
design  purposes  where  such  stage  parameters  may  not  be  available, 
Lieblein  [Ref.  11]  has  derived  an  equivalent  diffusion  ratio,  D  ,  which 
is  based  on  inlet  and  outlet  conditions  to  the  compressor  and  the  type 
of  blading  installed.  Unlike  a  single  airfoil  where  angle-of-attack  is 
used  to  predict  stall,  the  critical  variable  for  a  cascade  of  blades  is 
the  flow  turning-angle.  The  exit-flow  angle  as  well  as  the  inlet  angle 
must  now  be  considered. 

2.  Secondary  Flows  in  Axial  Compressors 

Returning  now  to  normal  compressor  operation,  flow  patterns 
based  on  axisymmetric  flow  distribution  and  an  infinite  number  of  blades 
have  been  described  quite  well  in  the  literature.  However,  it  is  necessary 
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to  consider  a  finite  number  of  blades  since  there  are  wall  boundary- 
layer  interactions,  tip  clearance  leakages,  blade  boundary-layers,  and 
centrifugal  effects  which  are  not  represented  adequately  by  assuming 
axisymmetric  flow.  Although  these  asymmetric  phenomena  have  been 
studied  independently  for  reasons  of  simplicity,  it  must  be  remembered 
that  their  effects  are  not  necessarily  additive  (non-linear)  and  may 
have  strong  interactions. 

Figure  8  illustrates  the  various  secondary  flows  which  have 
been  found  to  exist  in  axial-flow  turbomachines  as  a  result  of  the  above 
mentioned  viscous  and  three-dimensional  effects.  Lakshminarayana  and 
Horlock  [Ref.  12]  have  made  a  thorough  review  of  secondary  flows  in  axial - 
flow  turbomachines  and  have  indicated  where  in  the  literature  more 
specific  information  can  be  found.  Woods  [Ref.  13]  presents  the  current 
state  of  the  art  in  analytical  treatment  of  secondary  flows  and  associ- 
ated losses  in  axial-flow  turbomachines. 

The  induced  velocities  caused  by  secondary  flows  are  considered 
losses  since  only  a  small  part  of  their  kinetic  energy  is  recoverable 
downstream.  Secondary  velocities  may  perturb  the  exit  flow  angle  from 
a  blade  row.  Near  annulus  walls  where  lies  the  origin  of  the  central 
vortex  in  Fig.  8,  the  induced  exit  flow  angle  can  cause  local  blade 
stall.  Much  work  has  been  performed  in  analyzing  these  phenomena,  and 
some  analytical  methods  for  predicting  three-dimensional  exit  flow 
angles  have  been  proved  valid  through  experimental  results  [Refs.  14 
and  15]. 

Another  source  of  strong  secondary  flows  in  blade  passages  lies 
in  non-uniform  inlet  flow  such  as  described  previously.  Treating  inlet 
flow  distortion  as  vorticity  provides  a  natural  basis  from  which  to 
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FIG.  8        SECONDARY   FLOW  AND  VORTICES     IN  AN    AXIAL 
FLOW     COMPRESSOR      ROTOR     [adaptedfromRef.il] 

compute  resultant  secondary  circulation  at  the  exit  plane.  The  effects 
of  this  secondary  flow  on  cascade  performance  can  be  studied  analogously 
to  boundary-layer  induced  circulation.  Comparison  with  an  appropriate 
blade-loading  criterion  may  predict  the  compressor's  unstable  operation. 
3.  The  Theory  of  Vorticity  Analysis 

The  thermodynamic  properties  which  are  significant  in  compressor 
performance  are  total-pressure  (velocity),  total -temperature  (enthalpy), 
specific  heat,  and  molecular  weight.  Variation  of  any  one  or  combination 
of  these  properties  with  time  and/or  position  at  the  compressor  face 


25 


constitutes  inlet  flow  distortion  and  should  be  an  integral  part  of  any 
universal  distortion  parameter.  To  date  distortion  factors  dealing  only 
with  one  inlet  property  at  a  time  have  been  reported,  i.e.  total  - 
temperature  variation  ->  K-r  ->  Kj(stall)  or  total-pressure  fluctuation  ■*■ 
Kp  ■*■  Kp(stall).  This  type  of  piecemeal  approach,  while  providing  satis- 
factory results  for  certain  restricted  applications,  cannot  be  expected 
to  remain  valid  under  all  conditions.  In  reality,  an  aircraft  engine 
will  encounter  all  of  these  fluid  distortion  properties  at  one  time  or 
another.  A  partial  list  of  causes  of  these  is:  inlet  characteristics, 
environmental  conditions,  rocket  exhaust  from  ordnance,  hot  runway  opera- 
tion, steam  catapult  launches,  and  the  normal  operating  envelope  of  the 
aircraft.  There  is  definitely  need  for  a  more  realistic  or  universal 
distortion  parameter,  one  which  accounts  for  the  total  inlet  distortion. 
It  is  the  intent  of  this  work  to  show  that  a  distortion  parameter  based 
on  vorticity  would  be  more  flexible  and  yield  more  comprehensive  surge 
prediction  criteria  than  present  methods. 

Crocco's  Theorem  provides  a  versatile  means  of  converting  exist- 
ing inlet  distortion  maps  to  vorticity  in  that  pressure,  temperature, 
and  time  are  all  available  as  desired  for  inclusion  in  the  transformation. 
As  indicated  previously,  time  is  an  important  parameter  in  such  an  analy- 
sis, the  aspects  of  steady-state  versus  turbulence  induced  stalls  being 
quite  distinct.  Hawthorne  [Refs.  16  and  17],  Preston  [Ref.  18],  and 
Horlock  [Ref.  15]  all  deal  with  the  passage  of  vorticity  through  a  cas- 
cade. Horlock  and  Hawthorne  calculate  induced  exit  flow  angles  which 
involves  a  solution  to  Poisson's  Equation.  Calculating  appropriate  exit 
circulation  and  comparing  it  with  blade- lift  circulation  is  a  far  simpler 
process  that  yields  a  qualitative  feel  for  the  relative  size  of  this 
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effect.  Comparison  of  this  secondary  circulation  with  that  caused  by 
three-dimensional  effects  should  also  be  made  in  order  to  look  for 
possible  interactions. 

The  effect  of  inlet  distortion  on  blade  angle-of-attack  must 
also  be  considered  since  this  is  the  other  "half"  of  the  flow  turning 
angle.  This  effect  is  much  less  subtle  and  can  be  dealt  with  in  a  more 
straight-forward  manner.  Areas  of  low  pressure  represent  lower  axial 
velocity,  hence  higher  angles-of-attack.  The  unsteady  response  of  the 
blades  must  be  considered,  particularly  for  circumferential  inlet 
distortion  where  blade  passing  frequency  is  a  factor. 
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III.  CONVERSION  OF  PRESSURE  TO  VORTICITY  MAPS 


A.  CROCCO'S  THEOREM 

In  general  form,  Crocco's  Theorem  can  be  written  as 


T  grad  S  +  u  x  curl  u  =  grad  H  +  ^jf- 

9t 


(1) 


which  holds  in  unsteady  flow  [Ref.  19].  From  basic  Thermodynamics 
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Combine  Eqns.  (1)  and  (2),  assuming  H  =  Constant 


-5~  grad  PT  +  u  x  curl  u  =  -°— 


(3) 


The  following  development  will  apply  to  the  system  defined  in  Fig.  9, 
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Writing  Eqn.  (3)  in  component  form  using  cylindrical  coordinates  yields 


-RT  8PT  .  3ur  A 

e  +  (u„a)  -  iua)A)e   =  e  (4) 

p   3   r    e  z   z  eycr     .   r  v  ' 


•RT  1  3?T  „  3u 


e„+  (u„«_  -  u„a)7)eQ  =  e  (5) 


PT  r  39   D    z  r    r  Z  v     3t 


•RT  aPT  „  3uz  „ 

-ez  +  Kwe  "  Vr)ez  =  — '  ez  (6) 


PT  3z  fc  3t 


The  following  reference  quantities  taken  at  the  compressor  face  will 
be  used  to  non-dimensional ize  the  equations.  Refer  to  Fig.  9. 

Pj  =  Average  Stagnation  Pressure 
a  =  Average  Speed  of  Sound 
D  =  Compressor  Diameter 

w  =  Fluid  Vorticity 

Table  I  lists  some  of  the  characteristic  times  available.  Blade  passing 
frequency  will  be  used  in  this  development,  since  most  available  data 
include  compressor  speed;  then 
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The  numerical  value  of  the  reduced  frequency  for  stall  is  based  on  the 
work  of  Carta  iRef.  20].  For  values  less  than  approximately  0.1,  blade 
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TABLE       I                CHARACTERISTIC       TIMES 

BLADE   PASSING    FREQUENCY 

VT|pN          (1500)  (36) 

2ttR             (6)  (2) 

-  tuww  nz 

ONCE-PER-REV 

VT,P        _       ,50°     . 

125  Hz 

2ttR             (6)(2) 

WAVE    PROPAGATION 

RADIAL           -£L  s    ,20°  . 
2R          4 

300  Hz 

LONGITUDINAL               °*U 

=     60,180        =     120   Hz 

STALL    TIME 

REDUCED               u         wfc 
FREQUENCY                        U 

-  0.1 

f 

4(0.1)0200) 

(3)  (0.2)          200  Hz 
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stall  is  essentially  equivalent  to  steady-state  stall.  For  a  reduced 

frequency  larger  than  0.1,  stall  is  suppressed  due  to  dynamic  effects. 

Other  work  in  this  area  has  been  performed  by  Carta  [Ref.  22].  The 

t 
paper  by  Ericsson  and  Reding  [Ref.  21]  is  a  survey  of  dynamic  stall. 

Rewriting  the  component  equations  in  non-dimensional  form  and 

canceling  terms  yields 

1  aP'   a2  VTTDNa  /3u  \ 

yP'  3r'   a2     z    z       27ia2   \3f  / 


a2  VTIpNa  /3uA' 


+  -=■  (u   'a,  '   -  u  '«  ')     =    ~^r- i-^-]  (8) 


yPV   36        a2      Z     r  r    z  27ia2 


f) 


1     3P'       a2  VTTDNa,3u'x 

-(u/%'  •  uQ'w/)  =      -(-f-)  (9) 


YP'  32'   a2   r  e     6  r      2,a2 


Equations  (7),  (8),  and  (9)  are  exact.  Orders  of  magnitude  of  various 
terms  will  be  examined  to  ascertain  whether  or  not  complexity  can  be 
reduced  in  a  logical  and  consistent  manner. 

B.  ORDER  OF  MAGNITUDE  ANALYSIS 

Figure  A-2  represents  an  inlet-flow  distortion  map  at  the  face  of 
a  J85-GE-13  turbojet.  This  distortion  map  triggered  surge.  The  map 
of  Fig.  A-3  is  a  similar  situation  in  a  TF-30  turbofan  engine.  These 


o 

Reference  to  figures  and  equations  labeled  with  a  letter/number 

combination  have  more  detailed  explanations  which  are  contained  in  the 
Appendix. 
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maps  will  be  considered  representative  of  stall-inducing  flow  distortion 
for  purposes  of  making  an  order-of-magnitude  analysis  on  Eqns.  (7),  (8), 
and  (9). 

Establish  an  order  of  magnitude,  e  ,  such  that 

FT  "  PT 

L  .,  1  -  P'   e  0(e) 


PT 


From  Figs.  A- 2  and  A- 3 
-0.2  <  e  <  0.2 


Between  12  and  5  o'clock,  Fig.  A-2,  8P'/3r'  is  a  maximum.  Refer  to 
Appendix  A  for  all  numerical  values. 


3P'     3(P  /P  )     1.17  -  .97 

=  1.53 


3r'     3(2r/D      (7.1  -  6.05)(2/16.1 ) 

From  Fig.  A-3  at  the  12:30  o'clock  position 

3P'  1.04  -   .96 

.685 


3r'  (2.15  -  1.8)(2/6) 

3P' 

-0.7  <    <_  1.5 

3r' 


Between  5  and  7  o'clock,  Fig.  A-2,  a  maximum  value  of  3P'/3e  exists  such 

.416  per  radian 


that    3P'     1-20  -  .95 


39       .6  rad 
Similarly  from  Fig.  A-3  at  1  o'clock 
3P'     1.04  -  .92 


39       .35  rad 

■£1   ..4 


.343  per  radian 
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A  discussion  of  the  term  SP'/9z'  is  contained  in  Appendix  A,  Part  II 
An  acoustic  analysis  taken  at  the  compressor  face  yields  Eqn.  (A-2) 
which  allows  a  magnitude  estimation  of  aP'/az'. 
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From  Figs.  A-3  through  A-8  at  the  point  "0",  the  following  data  are 
available. 
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An  average  value  of  u  '  for  the  J85  can  be  obtained  from  data  in 

Appendix  C. 

u     460 

u  '  =  -L    =  =  .415 

z     a    1110 


A  maximum  value  would  be 

1/2 

lPT  "  Ps) 
P'  =  1.17       u7  *   VJ ±'  =  750  ft/sec 

p/2 

750 

u  '  - =  .675 

1110 

uz'  =  .7 


Data  with  which  to  estimate  the  magnitudes  of  ur'  and  uQ'  or  their 
partial  derivatives  are  not  available.  However,  an  upper  bound  of  order 
e  can  be  definitely  established  by  considering  the  consequences  of  ur' 
or  u  '  ~  0(e°).  A  radial  velocity  component  of  order  unity  is  not 
possible  within  a  duct.  Any  u  '  components  present  at  the  compressor 
face  are  caused  solely  by  turbulent  fluctuations  and  are  probably  of 
order  e  ,  however  an  upper  bound  of  e  will  be  assumed.  A  velocity 
component  such  as  u  ' ,  called  pre-rotation,  is  often  generated  intention- 
ally as  a  means  of  establishing  a  desired  angle-of-attack  on  the  first 
blade  row.  The  magnitude  of  this  velocity  should  be  smaller  than  u  ' 
and  certainly  not  larger.  A  velocity  in  the  circumferential  direction 
of  order  e  ,  about  equal  to  the  axial  velocity,  is  unrealistic;  and  if 
it  were  present,  the  engine  would  just  not  run.  For  the  case  under 
consideration  there  is  no  pre-rotation  so  that  any  u0'  which  does  exist 
is  caused  also  by  turbulent  fluctuations  and  is  probably  of  order  i   . 
Again  an  upper  bound  of  e   will  be  assumed. 


34 


Therefore 


Ufl'  a      u  '   <  U  ' 

o      r      z 

Partial  derivatives  of  u  '  and  u„'  with  respect  to  time  can  also  be 

o       r 

given  only  a  bound  with  respect  to  u  ' . 
au0'     3ur'     du   ' 


dV  dV  dV 

There  are  two  phenomena  causing  a  time  dependency.  One  is  due  to  the 
motion  of  the  rotor.  For  patterns  of  large  extent  on  the  Pj  map,  the 
appropriate  frequency  is  once-per-rev.  For  patterns  of  limited  extent, 
blade  passing  frequency  is  appropriate.  Unsteady  aspects  of  the  flow 
cause  a  time  dependency  also. 

Continuing,  the  magnitude  of  wr'  will  be  assessed. 


i 


ayD     /  8U       3uA  D 
u  '  e  __  =  (  —i.  + 

2a    \  ae    /dz  I  2a 


From  Fig.  A-2  between  8  and  9  o'clock,  and  data  from  Appendix  C 

(?  -  p  V/2 

P'  =  .97       u7  =  \— V  =  400  ft/sec 

p/2 


P'  =  1.17      uz  =  750  ft/sec 


A6  =  30°  =  .52  radians 

(350  ft/sec) (1.34  ft) 

wr'  = ~ = . =  .814  per  radian 

(.52  rad)(.605  ft)2(1110  ft/sec) 

wr'  -  .8 
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Similarly,  approximating  u^'  yields 
w0D    /  3un/  3u,  \  D 


/  3U^/  3UZ  \ 

\  p.  3r  / 


We      2a    V  ft.  ar  /  2a 


From  Fig.  A-2  between  12  and  5  o'clock 

-(750  -  400)  ft/sec  \   1.34  ft 


/-(750  -  400)  ft/sec  \ 
V      .06  ft      / 


=  -3.52 


.06  ft      /  2(1110  ft/sec) 


V  H  3'5 


For  the  axial  component  of  vorticity 


wzD    /  3u_   3u, 


wz 


\  3r    36  / 


2a    \  3r    36  /   2a 

J 


Since  both  partial  derivatives  are  really  unkown  in  size,  this  term 
will  be  evaluated  later. 

The  magnitude  of  3uz'/3t'  in  Eqn.  (9)  can  be  estimated  by  again 

referring  to  the  point  "B"  in  Figs.  A-3  through  A-8. 

1/2 
[  PT  -  Pc  )     '/*         /   (.92)(1610  -   1440)  \ 
.92 


(PT  "  Ps)    1/2         /  (.92)0610  -  1440)\ 
z  p/2  \  .0015/2  / 


uz     =     231   ft/sec 


P'     =     .96  uz     =     374  ft/sec 

At     =     (6  figs.)  x  (.2xl0~3  sec/fig.) 

(Auz\  3.14  /  j/4  -  C6\  v 

At  /  (1190)066)   \  1.2xl0"3     / 


3uz'  tt        /Au,\  3.14  /  374  -  231 


3t'  a  RPS  N 


3uz'   _     .06 


TV 
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Entropy  can  be  expressed  as 


S  -  S 
R 


=  -  Ln 


l_rj  y  - 


Ln 


B 


where  P  and  T  are  static  values  at  the  compressor  face.  Assuming  P  =  P, 
combination  of  this  equation  with  Eqn.  (2)  yields 


(+r  ■ + ■ 


Y-l 

(P'PT 


-   i+*r-e 


Y-l 
1  +  (e)  Y 


1  + 


The  deviation  of  local  speed  of  sound  from  average  speed  of  sound  is  less 
than  the  deviation  of  local  stagnation  pressure  from  average  stagnation 
pressure.  Therefore  a/a  will  be  set  equal  to  1. 

Table  II  is  a  summary  of  these  established  relative  magnitudes  in 
terms  of  E.  The  velocities,  ue'  and  ur' ,  are  placed  at  their  absolute 
upper  bounds  and  uz'  is  set  equal  to  e  for  later  verification. 


TABLE      H 


RELATIVE    ORDERS     OF     MAGNITUDE 


e2 

€ 

0 

6" 

P' 

dp 

dP' 

aP' 

it 

68 

JT 

a  -a 
a 

i         i 
r        e 

i 
uz 

OJr' 

OJg          OJZ' 

dur'    au; 

auz 

at'    at' 

at' 
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Now  Eqn.  (7)  can  be  written 


1 

yP' 

3P' 

ar' 

+  (ueV 

-     uz'V)  = 

VTIp  N  aur' 
2n     a  8t' 

c° 

e° 

c     e° 

e°c° 

£°          e2 

(7a) 


Dropping  terms  of  order  e 

1    3P' 

-  U-'wq'  =  0 

yP'  ar'  .     (10) 

-1   /3P' 


9      yP 


— (-) 


Similarly  for  Eqn.    (8) 

1         3P'  VTIp  N  3u   ' 

+     (uz'u>r'     -     ur'Uz')     =    —--^r  (8a) 


vP'r'     36  r  2tt    a  at' 

0  0  0     0  0  0         2 

ee  ee  ee  ee 


ay 


'u/r'Vae  / 


(ID 


The  size  of  oo   '  was  assumed  to  be  of  order  e     although  at  this  point  it 
is,  in  fact,  unknown.     In  arriving  at  Eqns.    (10)  and   (11),  wz'   is 
multiplied  by  small   terms   (u   '   and  u   ' ) ;  therefore  little  if  any  error 
should  be  introduced  here.     Subsequent  analysis  establishes  the  size  of 
wz' .     From  Eqn.    (9) 

I  9P'         /  \  V         N  9U    ' 

+(ur'V     -     ueV)   =    — -  —  (12) 

yP'     3z'        V        e  r  /  2tt    a  at' 

II  II  II  II 

0  0  0  0 

E  e  £        E  £        £  £  E 
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Equation  (12)  contains  most  of  the  terms  whose  magnitudes  were  generally 
uncertain.  All  terms  are  of  the  order  e  so  no  consistent  simplification 
can  be  made.  This  equation  will  not  be  used  in  the  subsequent  trans- 
formations from  Pj  to  vorticity. 

C.  APPLICATION  OF  RESULTS 

Equations  (10)  and  (11)  are  the  required  relations  for  transforming 
total -pressure  maps  to  circumferential  and  radial  vorticity  maps 
respectively.  The  procedure,  as  applied  to  the  J85  map,  is  outlined 
below. 

The  expressions  3P'/9r;  in  Eqn.  (10)  and  dP'/dd   in  Eqn.  (11)  are 
obtained  in  a  similar  fashion  from  Fig.  A-2.  For  one  particular  radial 
line,  values  of  P'  are  plotted  against  r',  the  slope  at  any  point  being 
the  local  value  of  3P'/3r'.  This  is  done  for  sufficient  radials  to 
adequately  cover  the  entire  map.  Local  values  of  P'  can  be  taken  direct- 
ly from  the  Pj  map,  and  from  these  u  '  can  be  calculated  (assuming 
constant  static  pressure).  Similarly  in  Eqn.  (11),  for  one  discrete 
radius,  values  of  P1  at  a  position,  9,  are  taken  at  intervals  through 
360  degrees.  These  values  are  plotted,  the  local  slopes  being  numerical- 
ly equal  to  SP'/9^.  Other  values  necessary  to  compute  wp'  can  be  taken 
directly  from  the  map  or  readily  computed.  See  Appendix  B  for  the 
calculations  outlined  above. 

The  results  of  these  calculations,  depicted  in  Figs.  10  and  11  as 
o)r'  and  w  '  contour  maps,  present  a  two-dimensional  picture  of  conditions 
at  the  compressor  face.  Since  v-(vxw)  =  0  ,  it  follows  from  Gauss' 
Divergence  Theorem  that 

->   ->  ->   -> 

ff   (oj  .  r)  dA  =  ///  (v  •  w)  dV  =  0 
A  y 
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Application  of  this  theorem  to  the  maps  of  u  '  and  u0'  contours  yields 
a  quantitative  estimate  of  the  axial  vorticity,  u>  '  which  is  reflected 
In  Table  II.  These  calculations  are  contained  also  in  Appendix  B,  and 
the  resulting  map  is  given  as  Fig.  12. 

Radial  and  circumferential  vorticity  maps  for  the  TF-30  total- 
pressure  data  in  Fig.  A-10  have  also  been  generated.  They  are  Figs. 
13  and  14. 
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Flow  Into  Paper 


Rotor 

Rotation 


x'dQ 


FIG.  10       RADIAL    VORTICITY    CONTOURS,        J 85 
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FIG.  II     CIRCUMFERENTIAL    VORTICITY    CONTOURS,  J85 
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FIG.    12 


CONTOURS     OF     NET     FLUX 
OF        Oj'z        FOR       J85 
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FIG.  13       RADIAL      VORTICITY      CONTOURS,      TF-30 
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FIG.  14 


CIRCUMFERENTIAL       VORTICITY 
CONTOURS,        TF  -  30 
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IV.  ANALYSIS  OF  VORTICITY  MAPS 

A.  RADIAL  VORTICITY  INLET  EFFECTS 

In  vector  notation  radial  vorticity  is  aligned  parallel  to  local 
blade  spans  of  the  compressor,  +cor  pointing  outward  radially  according 
to  the  right-hand  rule.  The  vorticity  in  Fig.  10  is  concentrated  in 
two  areas,  -a>r  between  5  and  7  o'clock  and  +cor  at  8  o'clock.  In 
analyzing  the  effects  of  these  local  concentrations  it  must  be  noted 
that  four  distinct  situations  can  be  represented  by  a  vorticity  map  such 
as  this.  These  are  listed  in  Table  III.  Vortex-sheets  are  aligned  par- 
allel with  the  flow  and  vortex  filaments  are  convected  with  the  flow. 

TABLE  III  Possible  Combinations  for  Radial  Inlet  Vorticity 


Case 

Vortex 

Blades 

Flow 

I 

Sheet 

Stator 

Steady 

II 

Sheet 

Rotor 

Unsteady 

III 

Filament 

Stator 

Unsteady 

IV 

Filament 

Rotor 

Unsteady 

Case  I 

A  continuous  sheet  of  vortex  filaments  passing  through  a  stationary 
blade  row  (stator)  is  a  steady-state  problem.  According  to  the  Biot- 
Savart  Law,  velocities  are  induced  in  the  fluid  which  cause  changes  in 
local  blade  incidences,  Fig.  15.  In  a  cascade  r  is  the  circulation 
about  a  blade.  I 


46 


and 


Then 


r    =     S<Vul   "  V 


L    =    pV  r 


AL  Ap  aV  Ar 

—    =    —    +    —    +    — 

L  p  Vm  r 


(13) 


(14) 


(15) 


r, 


LIFT 


A  Vaa 


FIG.   15         RADIAL    VORTEX  "SHEET     PASSING 
THROUGH      A      STATOR 


Therefore  blades  in  region  "A"  where  Ar/r  is  positive  are  more  heavily 
loaded. 

The  first  stator  row  is  behind  a  rotor  so  Figs.  10  and  13  cannot  be 
applied  here  without  knowing  the  influence  of  the  rotor  on  u'.     This 
case  could  be  applied  directly  to  inlet  guide-vanes  however. 
Case  II 

A  steady  vortex-sheet  passing  through  a  rotating  blade  row  is  an 
unsteady  problem  due  to  the  relative  motion  between  particle  paths  and 
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the  rotor.  Again,  induced  velocities  cause  variations  in  the  absolute 
velocity,  V-j ,  which  result  in  A3-j  (relative  inlet  flow  angle)  and  hence 
Ar  and  AL.  W  is  relative  velocity  in  Fig.  16.  As  blade  "A"  passes 
through  the  vortex-sheet  to  position  "B" ,  it  experiences  a  decrease  in 
lift  which  is  represented  by  a  smaller  riift-  This  results  in  shed 
circulation  equal  in  direction  to  rlif1.  and  a  resultant  decrease  in 
Strength  of  the  trailing  vortex-sheet.  Dynamic  response  of  the  blading 
is  a  factor  which  must  be  considered  in  this  case. 

For  a  rotating  cascade,  circulation  associated  with  blade  lift  is  a 
function  of  the  peripheral  velocity,  W  .  For  the  circuit  abed  in  Fig. 


16  this  is 


rlift  ■  S<Wul  -  Wu2> 
The  corresponding  lift  on  this  blade  is 

L  =  PrWro  =  PS(Wul  -  Wu2)Wo 


(16) 


(17) 


Shed  Vortices 

1 


ROTOR 


FIG.      16  RADIAL    VORTEX  "SHEET      PASSING 

THROUGH     A     ROTOR 


48 


For  our  sign  convention,  -u  causes  blade  unloading  while  +a>r  increases 
the  local  incidence.  Comparison  of  circulation  causing  blade  lift  to 
that  associated  with  inlet  distortion  (relist)  ^rom  ^9-  ^  should  give 
a  feel  for  the  relative  strength  of  this  effect.  For  brevity,  all  mean 
quantities  will  be  used  and  the  dynamic  response  of  the  blading  will  be 
neglected.  From  data  in  Appendices  C  and  D,  for  the  J85 

Wul  =  650  ft/sec 

Wu2  =  321  ft/sec 


1.095 


uu.  p 

»-m'      12 

Blade 

Chord  (C)  =  1,G5 
12 

riift 

= 

S(Wul  "  Wu2> 

Flift 

= 

30  ft2/sec 

rdist 

= 

^rcs 

2854.9 

Wp 

P' 

/9.17P'  -8.17 

ft 
ft 


(Relative) 


1       ,  3P' 


From  Fig.   A-2  at  8  o'clock 


P'  «  1 

avg 


From  Fig.   10  at  8  o'clock 

1_  /  aP'  \ 
r'  V  ae  J 


+  .13 


ae 

-1 


Therefore       ajr  =  +  372  sec        (absolute) 
and  rdist  =  2.97  ft2/sec    (absolute) 


(16) 


(18) 
CD-I) 
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The  ratio  of  circulations  becomes 


lift 


30 


10 


rdist     2.97 

which  is  about  10  per  cent.   It  must  be  remembered  that  this  is  occurring 
at  blade  passing  frequency,  and  any  induced  stall  would  be  dynamic  stall. 


D,  .  D   .   c         / 14,200  RevV31Blades\  ,  _ 
Blade  Passing  Frequency  =  ^~ 60~  Sec/V Rev/=  7'3 


Blades 
ms. 


Case  III 

Unsteady  vortex  filament  motion  through  a  stator  is  a  more  complex 
problem  since  time  dependency  of  the  fluid  motion  as  well  as  the  unsteady 
response  of  the  blading  must  be  considered.  This  situation  can  be 
pictured  as  the  passage  of  a  vortex  filament  or  a  local  concentration  of 
vorticity  passing  through  a  cascade.  In  a  similar  problem,  Crimi  and 
Reeves  [Ref.  25]  calculated  sizable  lift  fluctuations  as  a  helicopter 
rotor  blade  passed  over  a  spanwise  oriented  vortex.  Ward  [Ref.  26] 
observed  dynamic  stall  of  a  rotor  blade  which  he  asserted  was  caused  by 
a  radial  vortex  passing  beneath  the  blade.  Figure  17  depicts  this 
situation  in  a  cascade. 


FIG.    17  RADIAL    VORTEX    FILAMENTS 

PASSING    THROUGH    A    STATOR 
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Case  IV 

It  follows  that  a  vortex  filament  approaching  a  rotating  blade  row 
is  also  an  unsteady  problem.  Figure  18  shows  basically  how  Wu]  is 
affected  by  wr. 


-  OJ 


y, 


V 


Ql 


w2-n 


FIG.     18  RADIAL    VORTEX     FILAMENT 

APPROACHING       A      ROTOR 


B.  CIRCUMFERENTIAL  VORTICITY  INLET  EFFECTS 

Circumferential  vorticity  taken  vectorial ly  is  positive  in  the 
clockwise  direction.  Due  to  the  signs  in  Eqn.  10,  the  positive  gradients 
in  Figs.  11  and  14  represent  -|wq'|.  The  circumferential  vorticity  map 
for  the  J85  in  Fig.  11  shows  a  strong  concentration  of  vorticity  at  the 
outer  radii  extending  through  about  300  degrees.  Figure  14  shows  a 
similar  characteristic  for  the  TF-30  although  it  is  weaker  and  less 
extensive.  As  in  the  case  of  air',  four  basic  situations  can  exist  for 
possible  analysis.  These  are  listed  in  Table  IV. 
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TABLE  IV  Possible  Combinations  for  Circumferential  Inlet  Vorticity 


Case 


Vortex 


Blades 


Flow 


I 

Sheet 

Stator 

Steady 

II 

Sheet 

Rotor 

Steady 

III 

Filament 

Stator 

Unsteady 

IV 

Filament 

Rotor 

Unsteady 

Case  I 

Circumferential  vorticity  is  oriented  perpendicular  to  the  local 
blade  spans  of  the  compressor.  A  sheet  of  coQ  approaching  a  stationary 
blade  row  is  pictured  in  Fig.  19. 


Induced 
Velocity 


TIP 


FIG.     19 


CIRCUMFERENTIAL     VORTEX-  SHEET 
APPROACHING      A      STATOR 
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Induced  velocities  according  to  the  Law  of  Biot-Savart  cause  Va]  and 
hence  a-j  to  change  ahead  of  the  cascade.  This  results  in  variable  lift 
along  the  blade  spans  which  produces  shed  vorticity  downstream.  Shed 
vorticity  will  be  dealt  with  in  more  detail  later  when  considering 
secondary  flows. 
Case  II 

A  circumferential  vortex  sheet  approaching  a  rotating  blade  row 
presents  two  problems  for  analysis.  Relative  motion  between  the  particle 
paths  and  the  rotor  gives  rise  to  a  streamwise  component  of  vorticity  as 
well  as  a  component  perpendicular  to  the  blade  spans,  Fig.  20. 


FIG.  20     COMPONENTS  OF  CIRCUMFERENTIAL 
VORTEX- SHEET  APPROACHING  A  ROTOR 

The  concentration  of  w  '  at  the  outer  radius  in  Fig.  11  is  essentially 
constant  for  about  300  degrees  of  angular  displacement.  The  life  time 
of  one  particular  characteristic  such  as  this  can  be  as  long  as  the  time 
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for  one  revolution.  The  series  of  TF-30  PT  maps  in  Appendix  A  shows  a 
life-time  of  close  to  1/4  revolution.  The  significance  of  this  is  that 
one  blade  can  experience  the  same  velocity  distortion  for  up  to  360 
degrees  of  revolution  which  allows  a  steady-state  analysis.  Induced 
velocities  caused  by  the  normal  component  of  this  sheet  are  those  shown 
in  Fig.  19.  Velocities  induced  by  the  streamwise  component  are  pictured 
in  Fig.  21. 


Roior 


FIG.   21 


STREAMWISE      COMPONENT   OF    A 
CIRCUMFERENTIAL    VORTEX- SHEET 
AT     A     ROTOR 


The  streamwise  component  causes  the  typical  variations  in  lift  along  the 
blade  span,  but  it  also  passes  through  the  blade  row  and  affects  the 
exit  flow  angle.  This  will  be  discussed  later.  The  normal  component 
also  causes  a  spanwise  variation  in  lift  as  shown  in  Fig.  19. 

For  w0  of  small  circumferential  extent,  the  situation  is  unsteady 
due  to  relative  motion  between  particle  paths  and  the  rotor. 
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Case  III 

A  circumferential  vortex  filament  approaching  a  stator  is  an  un- 
steady problem.  The  effect  of  the  associated  entry  velocity  perturbations 
shown  in  Fig.  22,  a  function  of  the  dynamic  response  of  the  blading,  would 
probably  be  negligible.  However,  as  this  filament  passes  through  the 
blade  passage,  a  streamwise  component  develops  which  can  interact  with 
existing  secondary  flows.  This  phenomenon  will  be  discussed  in  the  next 
section. 


TIP 


duced  Velocities 


FIG.  22         CIRCUMFERENTIAL    VORTEX     FILAMENT 
APPROACHING        A        STATOR 


Case  IV 

A  circumferential  vortex  filament  approaching  a  rotor  is  different 
than  the  previous  case  only  because  a  streamwise  component  of  the  filament 
is  present  ahead  of  the  blades.  Again  unsteady  response  of  the  blades  is 
a  factor  to  consider. 
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C.  PASSAGE  AND  EXIT  EFFECTS 

Several  authors  have  treated  the  problem  of  turning  a  non-uniform 
flow  in  a  cascade  of  airfoils.  Squire  and  Winter's  method  for  calcu- 
lating secondary  flows  downstream  of  a  cascade  [Ref.  27]  has  been  used 
extensively  in  the  literature.  Preston  [Ref.  18],  using  a  more  simpli- 
fied approach,  also  treates  the  generation  of  streamwise  vorticity  which 
gives  rise  to  secondary  flow  after  a  cascade.  However,  Preston  assumes 
that  in  addition  to  the  components  of  vorticity  present  in  the  original 
vortex  filament,  a  change  in  circulation  about  the  blade  will  cause 
circulation  to  be  shed  into  the  stream.  Hawthorne  [Ref.  16  and  17] 
attempts  an  analytical  proof  of  the  existence  of  these  different  com- 
ponents of  vorticity  in  nonuniform  flow,  and  he  obtains  reasonable 
agreement  with  experimental  results. 

Consider  an  isolated  vortex  filament  in  a  sheet  normal  to  the  flow 
direction  as  it  passes  through  a  cascade  of  stationary  airfoils;  see 
Fig.  23.  Due  to  the  difference  in  flow  velocity  over  the  suction  and 
pressure  surfaces  of  the  blades,  this  filament  will  be  turned  and 
stretched  as  it  passes  through  the  cascade.  Downstream  of  the  cascade, 
segment  ab'  gives  rise  to  a  "distributed  streamwise  vorticity"  through- 
out the  passage  exit  plane.  This  description  applies  not  only  to  a 
filament  as  shown  here  but  also  to  a  vortex-sheet.  Squire  and  Winter's 
well  known  expression  for  this  distributed  vorticity  is 


*..  .  -  -2E£i  (19) 

dist       dz 


where  $  is  streamwise  vorticity,  e  is  the  deflection  angle,  and  Z  is 
the  spanwise  position. 
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Distributed 
Streamwise 
Vorticity 


FIG.   23 


Trailing  b 

Filament  Vorticity 


VORTEX      FILAMENT       PASSAGE 
THROUGH      A      CASCADE 


Segment  b-b',  part  of  the  original  filament  a-b  in  Fig.  23,  lies  along 
the  stagnation  streamline  and  in  considering  the  spanwise  direction,  this 
becomes  a  "vortex-sheet."  This  is  called  "trailing  filament"  vorticity. 
A  third  component  of  exit  vorticity,  "shed  circulation"  caused  by  vari- 
able spanwise  lift  due  to  the  asymmetric  inlet  flow,  also  lies  along  the 
stagnation  streamline  and  gives  strength  to  the  trailing  vortex-sheet. 
This  component  is  labeled  c-d  in  Fig.  23. 

Preston's  simple  analytical  expressions  for  these  secondary  components 
of  circulation,  as  applied  to  Fig.  23,  are  as  follows. 
Distributed  Secondary  Circulation: 

dr1  =  £2  s  cos  Cf)  dz 
And  replacing  ^   with  Squire  and  Winter's  expression  for  distributed 


secondary  vorticity,  Eqn.  (19)  ,  yields 

dFl  =  "2e(~dT")S  C0S  ^  dz 


(20) 
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Trailing  Filament  Circulation: 

dIV  =  i  dul\?  s  <;in  r__E_l 


•2  =&)2Ssi"<-r>  dz 


Trailing  Shed  Circulation: 


dr3 


.f^2S$1n  (-4-)  dz  (22) 

\  dz  /         c 


Trailing  shed  circulation  is  equal  in  magnitude  and  direction  to  dr2- 
The  net  circulation  associated  with  the  trailing  vorticity  in  one  blade 
passage  and  from  one  airfoil  in  the  form  of  a  "vortex  sheet"  is  the  sum 
of  these  three  components. 


dr  =  dr-|  +  dr2  +  dr3 
du 


dr 


=  (Jll)  [  4  S  sin  (|)   -  2  e   S  cos  (|)   ]  dz 


(23 


From  this  exit  secondary  circulation  the  change  in  mean  exit  flow  angle 
can  be  calculated.  This  total  circulation  is  not  uniformly  distributed 
over  the  exit  area  [Ref.  18]. 

Investigations  into  changes  in  deflection  due  to  cascade  secondary 
flows  have  generally  taken  two  lines  of  thought.  Values  of  induced 
deflection  at  the  lifting  line  were  derived  from  the  induced  velocities 
calculated  solely  from  the  strength  of  the  vortex  sheets  by  Carter 
[Ref.  28]  and  Hausmann  [Ref.  29]  to  mention  a  few.  These  methods  neglect 
the  effects  due  to  distributed  passage  vorticity  and  hence  are  unlikely 
to  provide  satisfactory  agreement  with  experiment  [Ref.  12]. 

Hawthorne  [Ref.  30],  Smith  [Ref.  31],  and  others  have  considered 
secondary  motions  to  be  induced  by  the  distributed  vortices  within  the 
passages.  It  is  assumed  that  the  flow  downstream  behaves  as  if  each 
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blade  passage  were  separated  from  its  neighbor  by  the  blade  wakes,  i.e. 
vortex  sheets.  The  secondary  velocities  in  the  blade  passages,  v 
circumferential ly  and  w  spanwise,  are  defined  by  a  stream  function,  \\>, 
such  that 

3^  3tf> 

W  =  -  ir~  V     =       7T- 

ay  9z 

The  stream  function  satisfied  continuity  for  this  steady  secondary  motion 
superposed  on  the  main  flow. 

iw  +  H     -  o 
3y    3z 

Substituting  the  stream  function  into  the  vorticity  equation   (w  =  curl   u) 
and  setting  it  equal   to  Squire  and  Winter's  expression  for  the  distributed 
vorticity,   Eqn.   19,  yields  Poisson's  Equation. 


l!*    +     *£*       =     .j       .     -Ze^l  (24) 

3y^  3z^  c  dz 


Solutions  to  Eqn.  (24)  are  tedious  but  have  been  accomplished  by  a  series 

method  [Refs.  27  and  30],  by  relaxation  techniques  [Ref.  17],  and  by 

electrical  analogs  [Ref.  32].  The  latter  reference  considers  variable 
spanwise  deflection  which  causes  Eqn.  (24)  to  become 

.    *2*  ■  -2e  &  +  ui  at  (25) 

An  exit  flow  solution  for  a  parabolic  inlet  profile  similar  to  the  case 
at  hand  is  given  by  Honda  [Ref.  33]. 
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Once  secondary  velocities,  v  and  w,  have  been  determined,  local 
changes  in  the  exit  flow  angle,  Aa2,  can  be  computed.  V2  is  the  main- 
stream exit  velocity  for  uniform  inlet  flow  with  exit  angle  o^. 


--U) 


Aap  =  Ta 


The  exit  flow  angle  for  distorted  inlet  flow  then  becomes 


dist 


Comparison  of  a2    with  appropriate  blade  loading  criteria  [Ref.  10 

dist 
and  11]  could  explain  the  compressor's  unstable  operation. 

Another  method  of  dealing  with  secondary  velocities,  v  and  w,  ob- 
tained from  solving  Poisson's  Equation  was  used  by  Hawthorne  [Ref.  17]. 
Again  the  assumption  of  blade  passages  extended  downstream  by  vortex- 
sheets  trailing  the  blades  is  made.  There  are  no  normal  velocities  at 
the  vortex-sheets.  The  strength  of  these  sheets  can  be  determined  by 
calculating  aw,  the  change  in  spanwise  secondary  velocity  across  them. 
As  mentioned  before,  these  sheets  of  fluid  leaving  the  trailing  edges 
are  vortex-sheets  composed  of  trailing  shed  and  trailing  filament 
vorticity;  see  Fig.  24.  If  the  trailing  filament  vorticity  is  estimated 
and  subtracted  from  the  total  strength  of  the  sheet,  the  trailing  shed  , 
circulation  remains.  From  this  can  be  calculated  the  variation  in  lift 
on  the  blade  since  shed  circulation  and  blade  lift  are  related  by 

AL  =  PUr3 

Hawthorne's  expression  for  this  trailing  filament  circulation  is 

j„            c/dul\M     sin  2a9  -  sin  2a-,  _  ,     /I        .  „, 

^     =     S{te)[z     COS  a,  •   ^   C°SM        (26) 
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FIG.  24       SECONDARY     FLOW    DOWNSTREAM     OF    A 
BLADE     ROW         adapted       from      Ref.    17] 


which  is  the  same  as  Preston's  formula  given  as  Eqn.  (21)  except  for 
different  inlet/exit  flow  angles.  Hawthorne's  predicted  shed  circulation 
agreed  quite  well  with  an  experimentally  determined  spanwise  lift 
variation. 

D.  CALCULATION  OF  SECONDARY  CIRCULATION 

A  calculation  of  these  components  of  secondary  circulation  as  caused 
by  the  circumferential  vorticity  in  Fig.  11  was  performed  in  order  to 
quantitatively  evaluate  the  effect  of  such  distortion  on  a  compressor. 
The  analysis  was  made  assuming  steady-state  conditions  which  means  Case 
II  of  Table  IV,  a  circumferential  vortex-sheet  passing  through  a  rotor. 
The  following  terms  require  definition. 

V  -  Absolute  Velocity 

W  -  Relative  Velocity 

U  -  Rotor  Peripheral  Velocity 
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a)  -  Absolute  Fluid  Vorticity 

C  -  Relative  Fluid  Vorticity 

fl  -  Rotor  Angular  Speed 

a  -  Absolute  Flow  Angle 

3  -  Relative  Flow  Angle 

•*■    -  Velocity  Vector 

-*-*-  -  Vorticity  Vector 

Subscripts;  see  Fig.  25 

r,e,z  -  Absolute  Coordinate  Direction 

s,n  -  Natural  or  Streamline  Coordinates 
s  -  Streamwise 
n  -  Normal  to  Streamlines 

0  -  Upstream  of  Compressor 

1  -  Compressor  face  (entry  to  blade  row) 

2  -  Blade-row  exit  plane 


Streamline 


FIG.  25 


NOMENCLATURE 
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Horlock's  method  [Ref.  15]  of  treating  the  passage  of  a  vortex 
filament  through  a  rotating  blade-row  results  in  the  following  vector 
transformation.  Figure  26(a)  is  the  velocity  diagram  ahead  of  the  rotor 
which  is  convecting  the  vortex  filament,  we  (one  of  many  in  a  sheet). 
The  absolute  velocity  is  axial  so  that  WQ-j  =  -urotor-  Figure  26(c)  is 
the  inlet  absolute  vorticity  combined  with  rotor  rotation  to  give  the 
inlet  relative  vorticity.  Equations  (27)  and  (28)  are  the  corresponding 
analytical  expressions  for  these  components  of  relative  vorticity  in 
terms  of  tuQ. 

Csi  =  +2  fi  cos  $i  -  ojQ  sin  3-j  (27) 

Cnl  =  -2  fi  sin  31  -  wQ  cos  B-|  (28) 

The  approach  of  a  vortex  filament  to  a  rotating  cascade  and  the  resultant 
secondary  relative  vorticity  at  the  passage  exit  are  pictured  in  Fig.  27. 
This  figure  is  similar  to  Fig.  23  except  it  applies  to  a  rotor  while  Fig. 
23  is  for  a  cascade. 

Combining  inlet  streamwise  vorticity  which  passes  through  the  cascade 
unaffected  with  the  exit  vorticity  of  Eqn.  (19)  yields  the  total  dis- 
tributed vorticity  in  the  blade  passage  exit  plane. 

^s2  =  Csl  "  2  e  Si  (29) 

Substituting  Eqns,  (27)  and  (28) 

CS2  =  (2fl  cos  3-|  -  w  sin  Bj)  -  2e(-wQ  cos  B-|  -  2ft  sin  b^)       (30) 

Since  this  exit  vorticity  is  in  the  direction  of  the  flow,  multiplying 
it  by  the  blade  passage  exit  area  normal  to  the  direction  of  flow  yields 
exit  distributed  circulation. 
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'2-1 


(c) 


FIG.    26         VECTOR      TRANSFORMATION       FROM 

ABSOLUTE     TO     RELATIVE     VORTIClTY 
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ROTOR 


FIG.  27         PASSAGE        OF       A       VORTEX    FILAMENT 
THROUGH    A    ROTATING     CASCADE 


drl 

dF    =    ^2  s  cos  62  (31) 

When  transformed  from  absolute  to  relative  coordinates,  Hawthorne's 
expression  for  trailing  filament  circulation  given  previously  as  Eqn. 
(26)  becomes 


di\ 
&r 


=     S 


^nl 


/  1  sin 


2B2  -  sin  23-j 


cos  e1 


-  2e  cos  B2 J   (32) 


Assuming  two-dimensional  flow,  Hawthorne's  equation  for  shed  circulation 
is 

dr0 


dr 


-S  cos  B-j  ^](tan  B2  -  tan  B-j) 


(33) 


A  simplified  application  of  these  equations  to  the  circumferential 
vorticity  and  total-pressure  maps  for  the  J85   (Figs.   11  and  A-2)  was 
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made  by  assuming  an  average  radial  distribution  of  aP'/Br"  and  P'  to 
exist  for  a  circumferential  segment  of  the  face.  Equation  D-2  was 
applied  to  radials  at  the  three  circumferential  positions  shown  in 
Fig.  28.  (Appendix  D) 


\p'/n 


-2854.9 


Vs 


7P'  -  8.17, 


ap1 

ar 


sec 


-1 


(D-2) 


Now  using  Eqns.  (27),  (28),  and  (29),  the  various  components  of  relative 
yorticity  were  calculated  for  a  sufficient  number  of  radial  positions; 
see  Appendix  E.  Figure  29  shows  the  relative  sizes  of  these  components 
for  circumferential  position  I. 


ap' 

dr' 


FIG.  28 


DATA    POINTS    FOR    CALCULATION    OF        tU 
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FIG.  29       VORTICITY    COMPONENTS    AT    RADIAL        I 
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Secondary  circulations  as  a  function  of  span  at  the  blade  row  exit  plane 
were  calculated  using  Eqns.  (31),  (32),  and  (33);  see  Appendix  E.  The 
area  under  curve  "A"  in  Fig.  30  represents  distributed  secondary  circu- 
lation. The  area  under  curve  "B"  is  trailing  filament  circulation,  and 
the  area  under  "C"  is  shed  circulation.  These  latter  two  are  the  cir- 
culations due  to  the  vortex-sheet.  The  magnitudes  of  these  components 
are  found  by  integration  over  dr'. 

.9 
R    /  (*E±\   dr'  =  32.1   ft2/sec 
.3  \dr  J 

where  R  is  the  radius  used  for  non-dimensional izing. 

Concerning  the  results  in  Table  V,  distributed  secondary  circulation, 
r-i ,  represents  a  loss  in  kinetic  energy.  Blade-lift  circulation  is 
approximately  equal  to  30  ft2/sec;  therefore  the  ratio  of  these  two  is 
about  1  to  1.  In  order  to  evaluate  the  actual  effect  r-j  has  on  the 
compressor,  a  solution  to  Poisson's  Equation,  as  described  previously, 
would  be  necessary.  From  this  solution  the  exact  strength  of  the 
trailing  vortex-sheet  could  be  determined.  Subtracting  trailing  fila- 
ment vorticity  (or  r2)  from  the  calculated  strength  of  this  vortex-sheet 
would  yield  the  exact  shed  circulation.  Equation  (33)  from  which  r3 
in  Table  V  has  been  calculated  was  found  by  Hawthorne  to  yield  values 
higher  than  those  found  experimentally.  However,  this  analysis  is 
essentially  an  order-of-magnitude  check,  and  Eqn.  (33)  is  adequate. 
Since  shed  circulation  is  related  to  the  spanwise  variation  in  blade 
lift  by 

AL  =  p  Wffi  r0 


then 


AL  =  rshed  (34) 


L    rlift 
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EXIT      SECONDARY       CIRCULATION 
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The  circulation  about  one  blade  in  the  rotor  under  consideration  for 
uniform  inlet  flow  is  approximately 


rlift  =  s<wul  "  Wu2>  =  30  ft2/sec 


The  average  value  of  positive  shed  circulation  over  the  three  circumfer- 
ential data  points  considered,  Radials  I,  II,  and  III,  is 


5.0  +  3.9  +  15.1 


=  8  ft2/sec 


So  the  variation  in  lift  from  Eqn.  (34)  is 

AL     8 

T  =  30   =   25  % 

which  represents  unloading  at  the  center  span  of  the  blade.  For  the 
negative  shed  circulation  at  the  blade  tip 


-6.0  -  8.5  -  10.9     AL     -8.8 

=  —  =  K  30  % 

3  L      30 


If  the  compressor  blading  is  operating  near  stall  it  certainly  seems 
plausible  that  a  change  of  25  to  30  per  cent  in  lift  could  trigger  stall 
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V.   SUMMARY  AND  REVIEW 

Figures  31  and  32  are  presented  as  a  review  and  outline  of  the 
preceding  theory  for  analyzing  the  effect  of  inlet  flow  distortion  on 
an  axial-flow  turbomachine.  Data  from  total -pressure  maps,  8Pj/3r  and 
3Pj/86,  consitute  the  initial  known  quantities.  All  pertinent  assump- 
tions and  simplifying  restrictions  are  included  in  the  chart  as  they 
apply.  The  resulting  determination  of  wz  through  a  consistent  order- 
of-magnitude  analysis  implies  complete  definition  of  the  velocity  field 
ahead  of  the  compressor  with  the  exception  of  the  ambiguity  in  sign  of 
w2.  This  is  significant  in  that  a  three-dimensional  flow  field  has  been 
determined  from  a  two-dimensional  pressure  map  at  the  expense  of  intro- 
ducing an  error  of  order  e^.  Continuing  in  Fig.  32,  the  analysis  next 
relates  portions  of  several  existing  theories  on  secondary  flows  in 
cascades  and  turbomachines  to  the  maps  of  radial  and  circumferential 
vortici.ty.  This  yields  an  estimate  of  the  strength  of  the  secondary 
flows  caused  hy  inlet  flow  distortion.  The  large  value  for  AL/L  which 
is  calculated  suggests  that  the  linear  cascade  theories  used  are  not 
adequate  for  the  accuracy  required  in  deriving  a  stall  criterion.  Some 
necessary  considerations  for  further  refinement  of  the  theory  which  were 
not  considered  in  this  analysis  are  (1)  blade  twist,  (2)  three-dimensional 
effects  and  interaction  between  normal  secondary  flows  with  distortion 
induced  flows,  and  (3)  changes  in  vorticity  with  passage  through  the 
blade  row. 
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Absolute  Vorticity,  w0  Map 


Relative  Vorticity  Components:  ^n-i ,  ^  -, 
Horlock  [Ref.  15] 


Ss2  ■  fUni) 

Squire  &  Winter  [Ref.  27] 


Total  Exit  Streamwise  Vorticity,  cS2 


Circulation  Components,  r 
Hawthorne  [Ref.  16] 
Preston  [Ref.  18] 


Li^  =  f^rshed) 
AL/L  =  Arshed/r 


Poisson's  Eqn.  ->  Induced  Velocity 
v  =  Normal  to  Span 
w  =  Spanwise 


w 


Vortex  Sheet  ->  aL 
Hawthorne  [Ref.  17] 


Change  Exit  Flow 
Angle,  A32 


Limiting  Blade  Loading 

Criteria,  Lieblein[Refs.   11   &  12] 


UNIVERSAL  DISTORTION  FACTOR 


Fig.  32.  Composite  Theory  for  Analyzing  Vorticity  Maps 
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VI.  CONCLUSIONS 

This  analysis  is  the  first  attempt  at  relating  total-pressure 
distortion  maps  to  the  internal  aerodynamics  of  a  turbomachine.  Sepa- 
rating the  inlet  distortion  into  components  of  radial  and  circumferential 
vorticity  introduced  flexibility  in  treating  the  problem  as  well  as 
maintained  a  comprehensible  physical  picture  of  the  aerodynamical 
situation.  For  the  data  analyzed  herein,  radial  vorticity  was  found 
to  be  less  significant  to  degradation  of  compressor  performance  than  was 
circumferential  vorticity  since  the  radial  vorticity  was  coupled  to  the 
dynamic  response  of  the  blading.  This  would  require  an  unsteady 
analysis.  The  relative  significance  of  u>r  and  w0  in  this  analysis  does 
not  necessarily  apply  in  general.  The  effect  of  axial  vorticity,  wz» 
must  be  included  in  further  work. 

The  size  of  the  calculated  shed  circulation  from  a  blade  indicates 
distortion  induced  secondary  flows  are  of  sufficient  size  to  seriously 
affect  the  compressor's  operation.  With  further  refinement  to  the 
present  theory,  secondary  flows  could  be  evaluated  to  a  degree  of 
accuracy  capable  of  predicting  stall  through  the  mechanism  of  shed 
circulation  or  changes  in  exit  flow  angles.  A  solution  to  Poisson's 
Equation  would  provide  the  needed  accuracy  with  which  to  use  appropriate 
blade-loading  criteria  in  defining  inlet-distortion  stall  criteria. 
Although  stall  normally  originates  somewhere  downstream  of  the  first  few 
stages,  the  distortion  must  be  treated  as  it  passes  through  each  indi- 
vidual stage.  The  present  theory  deals  only  with  the  first  blade  row; 
howeyer,  it  is  one  row  closer  to  the  origin  of  the  stall  and  has  elimi- 
nated a  certain  degree  of  uncertainty.  This  theory  may  be  extended  to 
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succeeding  stages  in  the  compressor  by  the  technique  illustrated  in 
Fig.  26. 

Finally,  in  working  toward  a  universal  distortion  factor,  the 
effects  of  other  types  of  inlet  distortion,  i.e.  total -temperature  and 
molecular  weight,  must  be  included.  These  extensions  to  the  present 
theory  should  provide  reliable  and  comprehensive  stall -inducing 
distortion  criteria. 
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APPENDIX  A 
ORDER  OF  MAGNITUDE  ANALYSIS 

Part  I  Preparation  of  Total -Pressure  Data 

The  following  data  and  calculations  apply  to  the  experiments  per- 
formed on  a  J85-GE-13  turbojet  as  described  in  Ref.  2,  a  =  0^  stall 
point.  Some  data  were  also  obtained  from  Ref.  24. 

Subscripts:  0  -  Free  stream 

1  -  Compressor  face 

2  -  Rotor  exit,  first  stage 
T  -  Stagnation  quantity 

FT1  =  1549  psf 

F|  =  1380  psf    (Static) 

TT1  =  532°R 

M  =  2.5 

00 

2 

Compressor  Inlet  flow  area,  A-j  =  185  in 

Figure  A-l  is  a  pressure  contour  map  of  Pji/Ptq  reproduced  from 
Ref.  2.  For  our  purposes,  this  must  be  converted  to  P'  which  is  defined 
as  PT1/PT1  . 

Pressure  Recovery  of  Inlet  =  .77 


FT1  . 

.77 

PT1 

_  PT1 

FT0 

PT1 

FT0 

FT1 

PT1 

PT1 
PT0 

1 
.77 

(A-l) 
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FIG.  A-l   PEAK  INSTANTANEOUS  PRESSURE 
CONTOURS,  J 85  (from  Ref.  2 ) 
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Application  of  Eqn.  (A-l)  to  Fig.  A-l  yields  the  desired  pressure 
contours,  Fig.  A-2.  Figures  A-3  through  A-ll  are  instantaneous  pressure 
contours  (P1)  taken  from  experiments  conducted  on  a  TF-30  engine  [Ref. 
23] .  The  radial  lines  at  the  centers  represent  compressor  rotation. 

RPM  -  10,000 

time  between  Figs.   =    /£0_      Sec\/JL      Rev  \ 

VlO4   Rev  A  32    Seg  / 

=  .1875  x  10"3 
«  .2  ms. 
•  Elapsed  time  (A-3  +  A-12)  «  1.5  ms. 
Figure  A-ll  is  the  frame  prior  to  engine  surge. 

Fs     =     1440  psf 

Estimated  Inlet  Values 
FT      =     .4 


pt  .  F    (wi^i  R2)^r 


FT    =     1610  psf 
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Flow  Into  Paper 


Rotor 
Rotation 


FIG.  A-2        TRANSFORMED        PRESSURE 

CONTOURS 
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FIG.  A-3 


FRAME     ONE     OF      INSTANTANEOUS 
PRESSURE     CONTOURS    SEQUENCE,    TF-30 
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FIG.     A-4 


FRAME       TWO       OF       INSTANTANEOUS 
PRESSURE      CONTOURS>     TF  -  30 
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FIG.     A- 5 


FRAME     THREE      OF      INSTANTANEOUS 
PRESSURE     CONTOURS    SEQUENCE,     TF-30 


83 


FIG.    A-6 


FRAME      FOUR      OF    INSTANTANEOUS 
PRESSURE    CONTOURS  SEQUENCE,  TF-30 
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■ 


FIG.    A-7        FRAME     FIVE     OF    INSTANTANEOUS 

PRESSURE    CONTOURS    SEQUENCE,  TF  -  30 
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FIG.   A-8 


FRAME    SIX    OF    INSTANTANEOUS 
PRESSURE    CONTOURS    SEQUENCE,  TF-30 
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FIG.  A-9  FRAME    SEVEN      OF    INSTANTANEOUS 

PRESSURE    CONTOURS    SEQUENCE  ,  TF-  30 
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FIG.  A-IO 


FRAME     EIGHT     OF    INSTANTANEOUS 
PRESSURE    CONTOURS    SEQUENCE, TF-30 
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FIG.  A-ll  FRAME    NINE     OF    INSTANTANEOUS 

PRESSURE    CONTOURS    SEQUENCE,  TF -30 
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Part  II  Pressure  Fluctuations  in  the  Z  Direction 

In  the  case  of  steady  flow,  aPj  is  present  only  when  dissipation 
occurs.  Some  causes  of  such  dissipation  in  an  inlet  are  shock  waves, 
flow  separation,  obstructions,  and  boundary  layers.  These  examples 
are  illustrated  in  Fig.  A-12.  


Boundary 
Layer 


(a)   Oscillating  Shock  Wave 


Flow 
Separation 


<^c£$£*^- 


(b)      Wall    Effects 


C 


S 


y 


Boundar 


y  Layer 


Strut 


C 


/ 


(c)      Obstructions 


FIG.   A-12 


SOURCES     OF    DISSIPATION 
WITHIN         INLETS 
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The  stagnation  pressure  gradient  is  present  only  in  the  vicinity  of  dis- 
sipation. Therefore  in  steady  flow,  aP'/az'  is  zero  at  the  compressor 
face  if  the  dissipation  is  located  upstream  as  shown  in  Fig.  A-13. 


C 


Streom  Tube 


FIG.  A-13         PRESSURE     HISTORY     FOR 
STEADY        INLET      FLOW 
For  unsteady  flow,  pressure  gradients  are  again  caused  by  dissipation. 
If  there  were  a  standing  shock  wave  in  the  inlet,  unsteady  flow  could 
cause  this  to  oscillate  in  the  axial  direction.     This  is  called  "shock 
wave  twitter"  and  can  be  self  excited  in  oscillation.     This  is  depicted 
in  Fig.  A-14. 


Range  of 
Displacement 


Xo 

Max 


-c 


Shock  Wave 


FIG.   A-14 


SHOCK- WAVE      TWITTER 
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The  changing  geometry  seen  by  this  oscillating  shock-wave  causes  pres- 
sure pulses  to  be  generated  as  shown  in  Fig.  A-15.  A/A*  is  the  area 
ratio  as  a  function  of  x.  The  second  quadrant  is  local  Mach  number  as 
a  function  of  area  ratio,  A/A*.  The  third  quadrant  is  the  stagnation 
pressure  ratio  for  a  shock  at  M.  Finally,  the  fourth  quadrant  yields 
the  derivative  of  Pj  with  respect  to  x. 


M  ■* 


FIG.  A-15        TOTAL   PRESSURE    FLUCTUATIONS 
FROM   OSCILLATING     SHOCK-WAVES 


The  resulting  aPt/ax,  a  function  of  diffuser  design,  can  be  evaluated 
at  the  compressor  face  by  assuming  an  acoustic  analysis  with  no  wave 
reflections  at  the  compressor.  Figure  A-16  is  a  plot  of  the  pressure 
pulses  (waves)  in  the  x-t  plane.  The  slope  of  the  characteristics  is 
equal  to  the  reciprocal  wave  speed;  curvature  is  due  to  the  variable 
area  at  the  center-body. 
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(time) 


FIG.    A-16 


PRESSURE    WAVES     IN     THE     x-f    PLANE 
AS   A  RESULT  OF  SHOCK- WAVE    TWITTER 


Consider  the  three  points  labeled  1,  2,  and  3  in  Fig.  A-16 


aP'         P2'  -  p,' 


3Z 


•2     "   zl 


Since  P3'  =  P   ' 


if 


p3'  -pr 


p2-  -v 


at'      VV      v -*i' 


equating 


3P' 


at' 


Z2'    -    Z    '    \      3P' 


V  -  Y 


8z' 
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v  -v 
v  -  V 


(z2  -  z^U/D)     1 

(t3  -  t^Ua/D)    a 


(z2  -  z^ 


(a  +  u)  =  a'  +  u 


therefore 


3P' 


1 


3P 


az'    (a'  +  u')  V  at' 


and 


3  P. 


3z 


1 


(a  +  u) 


3P- 


3t 


(A-2) 


The  usefulness  of  Eqn.  (A-2)  can  be  demonstrated  by  application  to  a 

simple  inlet/engine  configuration  with  an  oscillating  shock  wave  present; 

see  Fig.  A-17.  For  brevity  the  following  assumptions  are  made 

(1 . )  a  =  Constant 

(2.)  Wave  impedance  of  compressor  matches  wave  impedance 
of  duct 

(3. )  U  =  constant 

(4.)     No  wave  dispersion 

L    A 


Shock 
Wave 


C 


Xo  +  L 


FIG.  A-17         NOMENCLATURE 
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The  shock  wave,  oscillating  within  +N,  produces  pressure  waves  which 
are  represented  as  a  sine  wave  function. 

x  -  Xq  =  N  sin  (wt) 

ax 


within  AN 


8t 


=   w  N  COS  (wt) 


9Pl 
3t 


ax 

at 


a  P. 


ax 


U) 


N  cos  (tot) 


(29. 


^0 


Using  a  wave  transit  time  equal  to  L/(a+u)  with  a  =  wave  speed  relative 
to  the  gas,  the  unsteady  pressure  at  the  compressor  face  becomes 


a  Pi 


at 


=  w  N  cos  w(t  - 


L  ) 

a+u' 


VL 


ax 


Combining  this  with  Eqn.  (A-2)  yields 


aP- 


Sz 


N  cos  oj(t  -  -=— )  /  aPT 
a+u  / 


face 


u  +  a 


ax 


Non-dimensional ized  the  preceding  equation  becomes 

»  ft'  ,    ]     -^ 

3P'     (h)/n)(flN/a)  cos  [  Q  {Z        a'  +  u'    a 


)]/3P'' 


az 


3x'i 


where  n  is  rotor  angular  speed. 

This  is  one  method  of  relating  pressure  fluctuations  at  the  compressor 

face  (aP'/3z')  with  the  pressure  waves  generated  by  inlet/shock-wave 

interaction. 
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APPENDIX  B 
CONVERSION  OF  PT  MAP  TO  VORTICITY  -  J85 

Pressure  data  taken  directly  from  Fig.  B-l  are  listed  in  Table  B-I 
From  these  data  the  curves  in  Fig.  B-2  are  generated.  Slopes  taken 
from  these  curves,  which  are  tabulated  in  Table  B-II,  provide  a  means 
for  constructing  the  circumferential  vorticity  map  in  Fig.  11.  A 
similar  procedure  is  performed  in  generating  the  radial  vorticity  map 
in  Fig.  10.  Data  are  contained  in  Tables  B-I I I  and  B-IV  and  Fig.  B-3. 
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300 


270' 


240 


-  90< 


r  Scale 


.3         .5         .7  .9 

■      '■'.',' 


4  .6         .8         1.0 


FIG.  B-l  PRESSURE     DATA    AND    SCALES    FOR   J85 
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Table  B-I     Pressure  Data  Taken  in  Radial  Direction  from  J85 

Total-Pressure  Map 


0° 


P' 


30° 
r'         P! 


60° 
r*         P' 


90° 
r'         P' 


120° 
r'         P' 


150° 


.4  .91 

.68  .97 

.72  1.04 

.78  1.10 

.81  1.17 

.92  1.17 

.95  1.04 


.47  .91 

.7  .97 

.74  1.04 

.79  1.10 

.81  1.17 

.92  1.17 

.95  1.04 


.46  .91 

.75  .97 

.79  1.04 

.82  1.10 

.84  1.17 

.92  1.17 

.95  1.04 


.44  .91 

.62  .91 

.72  .91 

.75  .97 

.79  1.04 

.82  1.10 

.85  1.17 

.92  1.17 

.95  1.04 


.71  .91 

.75  .97 

.79  1.04 

.82  1.10 

.84  1.17 

.92  1.17 

.95  1.04 


.47  .84 

.55  .84 

.68  .91 

.75  .97 

.78  1.04 

.80  1.10 

.82  1.17 

.92  1.17 

.95  1.04 


180° 
r'     P' 


210° 
r'         P1 


240< 


1         P ' 


270° 


P* 


300° 


P' 


330< 


.4  .84 

.6  .84 

.7  .91 

.78  .97 

.82  1.04 

.86  1.10 

.90  1.1Q 

.92  1.04 

.97  .97 


82 


91 


.42  .91 
.9  .97 
.95       .97 


.3  .91 

.43  .97 

.65  1.04 

.75  1.10 

.79  1.17 

.92  1.17 

.95  1.04 


.29  .91 

.47  .97 

.68  1.04 

.76  1.10 

.80  1.17 

.92  1.17 

.97  1.04 


.3  .91 

.6  .97 

.71  1.04 

.76  1.10 

.8  1.17 

.92  1.17 

.95  1.04 
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Table  B-II  Pressure  Gradients  Obtained  from  Slopes  in  Fig.  B-2 


0° 

30° 

60° 

90° 

120° 

150 

0 

r' 

3P7ar' 

r'   3P'/3r' 

r'  3 

P'/3r' 

r'   3P'/3r' 

r'   3P73r' 

r'   3P'/3r' 

.4 

0 

.4 

.13 

.4 

.07 

.32 

.045 

.4         .03 

.3 

-.22 

.5 

.09 

.5 

.16 

.5 

.12 

.38 

.17 

.5         .06 

.3 

-.22 

.6 

.28 

.6 

.26 

.6 

.135 

.44 

.26 

.6         .1 

.5 

-.22 

.7 

.9 

.7 

.5 

.7 

.25 

.52 

0 

.68       .33 

.51 

0 

.76 

1.8 

.76 

1.4 

.74 

.56 

.6 

-.26 

.72     1.08 

.54 

.3 

.78 

2.6 

.79 

5 

.8 

2.3 

.62 

-.27 

.8         3 

.6 

.53 

.80 

1.6 

.8 

6 

.82 

4.5 

.66 

0 

.82     2.5 

.7 

.73 

.82 

.83 

.82 

1.1 

.86 

.55 

.7 

.29 

.84     1.08 

.76 

1.5 

.84 

.37 

.86 

0 

-.89 

0 

.72 

1.2 

.86       .4 

.79 

4.0 

.87 

0 

.9 

-.72 

.92 

-.6 

.75 

2.0 

.89       0 

.8 

7 

.9 

-.6 

.94 

-2.7 

.94 

-2.7 

.8 

2.0 

.9       -.26 

8.2 

1.3 

.94 

-2.7 

.96 

-3 

.96 

-3.0 

.83 

6.0 

.92     -.6 

.86 

0 

.96 

-3.0 

.85 

.87 

.9 
.2 

.94  -2.7 
.96  -3.0 

.9 
.94 

-.72 
-2.7 
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Table  B-II     (Continued) 


180 

0 

210° 

240' 

j 

270° 

300° 

330° 

r'  3P 

7ar' 

r'  3P'/3r' 

r'  3P 

73r' 

r'  3P'/3r' 

r'  3P'/3r' 

r'  3P'/3r' 

.34 

-.25 

.3    .07 

.3 

.13 

.32   .505 

.4    .35 

.4    .15 

.4 

-.35 

.3    .07 

.3 

.13 

.4    .38 

.5    .26 

.6    .33 

.48 

-.12 

1.0 

| 

.5    .28 

.6    .35 

.7    .88 

.51 

0 

.9 

.13 

.6    .35 

.7    .55 

.76  1.9 

.6 

.37 

.92 

0 

.7    .52 

.74   .8 

.78  2.1 

.7 

.7 

.94  • 

-.15 
i 

.76  1.6 

.78  2.3 

.8   1.08 

.78 

1.7 

1 

1 

.78   .56 

.8   1.3 

.82   .6 

1 

1 

1.0   - 

-.15 

.8    .44 

.82   .65 

.84   .23 

.86 

1.7 

.85   0 

.85   0 

.86   0 

.88 

0 

.89  -.37 

.9   -.7 

.9   -.72 

.9   • 

-2.5 

.92  -.6 

.92  -1.3 

.94  -2.7 

.96  • 

-1.2 

.94  -2.7 

.94  -2.7 

.96  -3. 

.96  -3.0 

.96  -3. 
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Table  B-III  Pressure  Data  Taken  in  Circumferential  Direction  from  J85 

Total -Pressure  Map 


r'  =  .38 

•  r'   = 

.46 

r'   = 

.54 

r'   = 

.61 

r'  = 

.69 

r'  = 

.77 

e         P' 
(rad) 

6 

(rad) 

P' 

0 

(rad) 

P» 

0 

(rad) 

P' 

0 

(rad) 

P' 

0 

(rad) 

P' 

4.26       .91 

4.19 

.91 

1.81 

.91 

1.69 

.91 

0 

.97 

0 

1.04 

6.17       .91 

4.48 

.97 

2.35 

.84 

2.41 

.84 

1.25 

.91 

.872 

.97 

5.23 

.97 

3.14 

.84 

3.19 

.84 

2.79 

.91 

2.09 

.97 

.174 

.91 

4.10 

.91 

3.98 

.91 

3.84 

.91 

2.96 

.97 

4.41 

.97 

4.36 

.97 

4.31 

.97 

3.40 

.91 

5.48 

.97 

5.72 

.97 

4.59 
5.13 

1.04 
1.04 

3.87 
4.27 
4.48 

.91 

.97 

1.04 

r'   = 

.85 

r'   = 

.93 

0 

(rad) 

P' 

0 
(rad) 

P' 

0 
(rad) 

P' 

.174 

1.17 

3.70 

v91 

3.03 

1.17 

.785 

1.10 

3.92 

.91 

3.19 

1.10 

1.39 

1.04 

4.22 

.97 

3.37 

1.04 

2.0 

1.10 

4.39 

1.04 

3.58 

.97 

2.26 

1.17 

4.52 

1.10 

4.10 

.97 

2.61 

1.17 

4.71 

1.17 

4.27 

1.04 

2.87 

1.10 

4.88 

1.17 

4.45 

1.10 

3.14 

1.04 

5.23 

1.17 

4.57 

1.17 

3.40 

.97 

6.11 

1.17 
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Application  of  Gauss'  Divergence  Theorem  is  performed  in  the  follow- 
ing manner.  A  grid  such  as  is  shown  in  Fig.  B-4  is  placed  over  a 
vorticity  map. 


FIG.    B-4 


METHOD     OF     DIVERGENCE 
APPLICATION 


The  net  a>  '  passing  through  region  abed  is  equal  to 


[(-.13)(*|)  +  (0=^11) (lb)]  -  [(^1|=0)(|  cd)  +  (0^11) (|  cd)]  = 


net 


Net  flux  of  wQ'  is  calculated  in  a  similar  manner  using  the  same  grid 
over  an  w0'  map.  The  sum  of  flux  of  »  '  and  w0'  is  then  equal  to 
w7'(ab-cd)  so  that 


w. 


to,'  (ab*cd) 


Area(abcd) 


This  procedure  is  performed  on  each  of  the  36  segments  so  that  a  complete 
o)Z'  map  can  be  drawn.  These  data  are  contained  in  Table  B-V.   If  the  sum 
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of  net  flux  of  u>r'  and  to'  is  negative,  then  the  net  flux  of  w  '  must 
be  positive.  There  is  an  ambiguity  of  algebraic  sign  which  cannot  be 
resolved  without  additional  data  or  information.  One  evaluates  w'-fi  , 
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APPENDIX  C 
ESTIMATION  OF  SOME  J85  PARAMETERS 

Subscripts:  T  -  Stagnation  Quantity 

0  -  Freestream  ahead  of  Compressor 

1  -  Compressor  Face 

2  -  Rotor  Exit,  first  stage 
m  -  Mean  radius 

Data  obtained  from  Ref.  24: 

FT1  =  1549  psf 

P-|  =  1380  psf 

TT1  =  532°R 
Data  obtained  from  Ref.  2  and  Ref.  34: 


518  7°R  (Corrected  Total -Temperature) 


a  k 


J (Corrected  Total-Pressure) 


2116  psf 
N*  =  16,500  RPM  (Rated  Engine  Speed) 


Wcorr   =  ^r-  =  15  KG/sec  =  33.5  lbs/sec    (Corrected  Mass  Flow  Rate) 


Percentage  Corrected  Engine  Speed  =  JL--/6  x  100  =  86.9 

D-j  =  16.1  inches  (Compressor  Face  Outside  Diameter) 

Dl  hub  =  4-8  inches 

A,  =  155  in2  (Entrance  Flow  Area) 

N  =  31  (Number  of  blades  in  first  row) 
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FfG.   C-l  BLADE       PASSAGE       DIMENSIONS 


Calculated  Parameters 


Tl 


JL 


(i  +rJ_M2)    Y-] 

Y 


1549 


3.5 


ffgj    -     1.125     =     (1   +.2MO 


M1     =     .414 


Tl 


=     1   +  3£l  M2 


Tl     =     TT1 


/U 


1m2 


T]   =     514°R 
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a  =    Y^ 


IL  =  M  a 


RT, 


a  =  49.02  V514 


a  =  1110  ft/sec 


IL.  =  (.414)(1110  ft/sec) 


460  ft/sec 


1380  lbs/ft2 


H715.6|^|)  (514°R) 


a 


=  .001565  slugs/ft3 


As  a  check  on  these  figures,  another  method  may  be  used 


1549 
w  =   WCOr  6  =  (33.5  lbs/sec)  j  2TT6  j 

yr         V532/518 


W  =  PA1UZ 


W  =  24.15  lbs/sec 


W 


pA- 


(24.15  lbs/sec)  (144  in2/ft2) 

(.001565  s1u9s  )  (155  in2)  (32.2  J^L   ) 
ft^>  slug 

U7  =  462  ft/sec 
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Estimation  of  3n 


W 


ui 


FIG.    C-2         VELOCITIES    AHEAD 

OF   ROTOR 


W  +  Relative  Velocity 
V  ■>  Absolute  Velocity 


&,     =     tan"YWul 


(M\    =     tan-Vw  -  D  •   RPMN 


DM     =     .875  ft 


Ang.   Speed     =     (86.9)(16500) 

VTT/518.7     x  100 


=     14,200  RPM 


Wnl     =    tt  (.875  _1L_  )   (Ji200     rev} 


rev 


60    sec' 


Wul  =  650  ft/sec 


Mean  Radius 


«   _  t;m-l  /  650  » 
31  "  tan   (  460  } 


3-]  =  55°  Mean  Radius 
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Estimation  of  32 

Individual  stage  characteristics  and  loss  data  for  the  J85  are  not 
available.  The  following  data  are  from  Ref.  34: 

Compressor  of  Free-Vortex  Design 

Lieblein  Diffusion  Parameter,  D  <  .5 

Solidity  (a)  =  1.14 


io  =  .96  J 


First  Stage  Rotor 
Axial  Velocity  Ratii 


D  =  (1  -  ^  )  +  ^i.  (Ref.  10) 

Wl     2VU} 

w2  =  va2/cos  32  =  ,96  Val/cos  32 
=  442/cos  32  =  546 

Wu2  =  Va2  tan  B2  =  .96  Vfll  tan  e2 
=  442  tan  62 

W-|  =  Va/cos  B-|   =  802  ft/sec 

5  =  (i    442/cos  62  \  +  (650  -  442  tan  32) 
802  (2)(1.14)(802) 

32  =  36°  minimum 

Wu2  =  320  ft/sec 

Since  the  loading  limit  is  D  <  .5,  36°  can  be  considered  a  minimum  for 
32>  and  the  maximum  turning  angle  becomes 
A3  =  3]  -  32  =  55°  -  36° 
=  19° 
e  =  .332   radians 
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APPENDIX  D 
CALCULATION  OF  ABSOLUTE  VORTICITY 


Radial  Vorticity: 

1      /  1   3P' 


r     yP'Uz'   V  r'  36  / 


(ID 


Y   =   1.4 

p,  m  !n 


PT1 


a     Uz/M       Uz 


m     R  ^(PT-PS)(2/P)\1/2  .  R  /PT-  1380  \1/2 

(Pt-PS)(2/p)/     v1549  -  l380y 


-  (   Pt         \  1/2 


1/2     ,     Pt      N  1/2 


=  pfffl.fl_  -  8.17)   "-  m(9.17  J-  -  8.17^ 
V  PT   169        '        V     PT      ' 


Uz'  =  M./9.17P'  -  8.17 

wr   =  — — 
2a 

Substituting  into  Eqn.  (11) 

B   -  2a"  /I  9P_|_  \  ft/sec 

r    1.4  D  P'  M  V  9.17P'-8.17  (  r"  86  j    ft 
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Recall:     a  =  1110  ft/sec 
D  -  Ifjl  ft 

M  =  .414 


03 


2854.9       /  i  3p'  v    _-, 


r     P'  -/9.17P'  -  8.16 


I1-—)   ^c-1  (D-D 


Circumferential  Vorticity 


---   -2854'9    -(^-\      sec"1  (D-2) 

1.16  \   ar-y 


P'  -/9.17P'  -  8. 


Equation  (D-2)  is  used  to  convert  the  contours  in  Fig.  11  to 
absolute  vorticity  for  subsequent  calculation  of  secondary  circulation. 
Tabulated  results  are  contained  in  Table  D-I. 
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Table  D-I  Absolute  Vorticity  for  Three  Circumferential  Positions 

on  »  '  Map 


Rad- 

ial 

I 

Radial 

II 

r'     / 

r  (in) 

P'         3P7: 

)r' 

u0   (sec-1) 

P' 

3P7sr' 

u>e   (sec-1) 

.3 

2.43 

.91            .: 

35 

-2554.9 

.9 

.05 

-520 

.325 

2.63 

.93 

-1770.8 

.05 

.35 

2.83 

.94 

-1567.6 

.05 

.375 

3.03 

.95 

-1416.3 

.05 

.4 

3.23 

.96 

-1297.8 

.05 

1 

1 

.425 

3.43 

.97 

-1201.6 

' 

■ 

.06 

-624.1 

.45 

3.63 

.98 

-1121.4 

.9 

.075 

-780.1 

.475 

3.83 

.99 

-1053.2 

.91 

.09 

-656.9 

.5 

4.03 

.99 

-1053.2 

.91 

.1 

-729.9 

.525 

4.23 

1.0 

-994.2 

.91 

.11 

-802.9 

.55 

4.43 

1.0 

-994.2 

.92 

.13 

-767.3 

.575 

4.63 

1.01 

-942.5 

.92 

.14 

-826.3 

.6 

4.83 

1.01 

-942.5 

.93 

.15 

-758.9 

.625 

5.03 

1.02 

-896.7 

.93 

.17 

-860.1 

.65 

5.23 

1.03              * 

i 

-855.7 

.94 

.20 

-895.7 

.675 

5.43 

1.04            .: 

35 

-818.8 

.94 

.23 

-1030.1 

.7 

5.63 

1.05         .; 

75 

-1682.7 

.95 

.25 

-1011.6 

.725 

5.83 

1.06           l.( 

) 

-2156.2 

.97 

.5 

-1716.6 

.75 

6.04 

1.08        i.; 

I 

-2402.2 

1.0 

1.0 

-2840.7 

.775 

6.24 

l.io        i.; 

I 

-2243.5 

1.04 

1.5 

-3509.2 

.8 

6.44 

1.17           l.( 

] 

-1522.4 

1.10 

2.0 

-3739.2 

.825 

6.64 

1.18         .; 

75 

-1112.4 

1.17 

1.5 

-2283.6 

.85 

6.84 

1.19            0 

0 

1.18 

.5 

-741.6 

.875 

7.04 

1.20        -.; 

>5 

352.6 

1.20 

0 

0 

.9 

7.24 

K18        -.; 

75 

1112.4 

1.18 

-.5 

+741.6 

.925 

7.44 

1.17         -l.( 

) 

1522.4 

1.17 

-1.0 

1522.4 

.95 

7.65 

1.14        -l.i 

2480.2 

1.1 

-3.0 

5608.9 

.975 

7.85 

1.1           -3.( 

3 

5608.9 

1.05 

-3.2 

7179.8 

1.0 

8.05 
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Table  D-I   (continued) 


Rad 

ial 

III 

p'  / 

r  (in) 

P' 

3P'/ar' 

ujq   (sec-1) 

.3 

2.43 

.9 

-.2 

+2080.3 

.325 

2.63 

.35 

2.83 

.375 

3.03 

.4 

3.23 

.425 

3.43 

.45 

3.63 

i 

.475 

3.83 

i 

■ 

1 

.5 

4.03 

-.2 

+2080.3 

.525 

4.23 

0 

0 

.55 

4.43 

+  .3 

-3120.5 

.575 

4.63 

.4 

-4160.7 

.6 

4.83 

.5 

-5200.8 

.625 

5.03 

i 

1 

.5 

-5200.8 

.65 

5.23 

.90 

.5 

.675 

5.43 

.90 

.5 

.7 

5.63 

.91 

.5 

-3649.9 

.725 

5.83 

.94 

1.0 

-4478.9 

.75 

6.04 

.97 

1.5 

-5149.8 

.775 

6.24 

1.0 

2.0 

-5681.4 

.8 

6.44 

1.04 

2.3 

-5380.8 

.825 

6.64 

1.17 

1.7 

-2588.0 

.85 

6.84 

1.18 

0 

0 

.875 

7.04 

1.2 

-.3 

+423.2 

.9 

7.24 

1.18 

-.75 

+1112.4 

.925 

7.44 

1.17 

-1.5 

+2283.6 

.95 

7.65 

1.04 

-3.0 

+7018.4 

.975 

7.85 

1.0 

-3.2 

+9090.3 

1.0 

8.05 

119 


APPENDIX  E 
•  CALCULATION  OF  CIRCULATION  FROM  u>6  MAP  FOR  J85 

C$i  =  +  2ft  cos  g,  -  w0  sin  3-j  (28) 

=  +  2  (236.6  j^)  cos  55°  -  u>e  sin  55° 

;S1  =  +  271.4  -  .8191  w0  (E-l) 

Cni  =  -  2ft  sin  3-|  +  o)e  cos  B1  (29) 

=  2  (236.6  M)  sin  55°  +  (cos  55°)  ufl 
Sec  D 

Cnl  =  -  387.6  -  .573  u»Q  (E-2) 

CS2  =  +  CS1  -  2e?nl  (29) 

=  (271.4  -  .8191  co0)  +  2(.3316  rad)(387.6  +  .573  u>Q) 
CS2  =  528.4  -  .439  Wfi  (E-3) 


=     ^52  s  cos  32  (31) 

dr 

s     =     1  095  in         m      m2b 
12  in/ft 

32  =  36° 

drl 

gp  =  .0738  ?S2    ft/sec  (E-4) 
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/ 1  sin  2B2  -  sin  26i    _     Q  \  ,__. 

=  ScJt  " L  "  2e  COS  32  (32) 

n  l\  c  rnc  R,  ' 


dr2 

dr~       'nl\ c      cos 


B1  =  55° 


e  =  .3316 


dFp 

^   =  .0480  ?nl   ft/sec  (E-5) 


dr 


S  cos  ^  C-.  (tan  g2  -  tan  S])  (33) 


dr, 

—L    =     .0336  r  (E-6) 

dr  ni 


Tabulated  results  for  Eqns.  (E-l)  through  E-6)  using  previously 
determined  data  are  contained  in  Tables  E-I,  E-II,  and  E-III. 
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Table  E-I  Relative  Vorticity  and  Circulation  Components  at  Radial  I 


r' 

O)0 

^Sl 

H 

CS2., 

diydr 

dr2/dr 

dr3/dr 

(sec-1) 

(sec~l ) 

(sec-' ) 

(sec-1) 

(ft/sec) 

(ft/sec) 

(ft/sec) 

.3 

-2554.9 

+2364.1 

+1076.3 

+1650 

+121.7 

51.2 

+36.3 

,325 

-1770.8 

+1721.8 

+627.0 

+1305.7 

+96.3 

29.8 

+21.1 

.35 

-1567.6 

+1555.4 

+510.6 

+1216.5 

+89.7 

24.5 

+17.2 

.375 

-1416.3 

+1431.4 

+423.9 

+1150.1 

+84.8 

20.1 

+14.3 

.4 

-1297.8 

+1334.4 

+356.0 

+1098.1 

+81.0 

16.9 

+12.0 

-1201.6 

+1255.6 

+300.9 

+1055.9 

+77.9 

14.3 

+10.1 

.45 

-1121.4 

+1189.9 

+254.9 

+1020.6 

+75.3 

12.1 

+8.5 

-1053.2 

+1134.0 

+215.8 

+990.7 

+73.1 

10.3 

+7.3 

.5 

-1053.2 

+1134.0 

+215.8 

+990.7 

+73.1 

10.3 

+7.3 

-994.2 

+1085.7 

+182.0 

+964.8 

+71.2 

8.6 

+6.14 

.55 

-994.2 

+1085,7 

+182.0 

+964.8 

+71.2 

8.6 

+6.14 

-942.5 

+1043.4 

+152.4 

+942.1 

+69.5 

7.2 

+5.12 

.6 

-942.5 

+1043.4 

+152.4 

+942.1 

+69.5 

7.2 

+5.12 

-896.7 

+1005.8 

+126.2 

+922.0 

+68.0 

6.0 

+4.2 

.65 

-855.7 

+972.3 

+102.7 

+904.0 

+66.7 

4.8 

+3.5 

-818.8 

+942.0 

+81.5 

+887.8 

+65.5 

3.9 

+2.75 

.7 

-1682.7 

+1649.6 

^+576.5 

+1267.1 

+93.5 

27.4 

+19.4 

-2156.2 

+2037.5 

+847.9 

+1474.9 

+108.8 

40.3 

+28.7 

.75 

-2402.2 

+2239.0 

+988.8 

+1582.9 

+116.8 

47.0 

+33.3 

-2243.5 

+2109.0 

+897.9 

+1513.2 

+111.6 

42.6 

+30.3 

.8 

-1522.4 

+1518.3 

+484.7 

+1196.7 

+88.3 

23.0 

+16.3 

-1112.4 

+1182.5 

+249.8 

+1016.7 

+75.0 

11.8 

+8.4 

.85 

0 

+271.4 

-387.6 

+528.4 

+38.9 

-18.4 

-13.0 

352.6 

-17.4 

-589.6 

+373.6 

+27.5 

-28.0 

-19.8 

.9 

1112.4 

-639.7 

-1025.0 

+40.0 

+2.9 

-48.8 

-34.5 

1522.4 

-975.5 

-1259.9 

-139.9 

-10.3 

-60.0 

-42.5 

.95 

2480.2 

-1760.1 

-1808.7 

-560.4 

-41.3 

-86.0 

-61.0 

5608.9 

-4322.8 

-3601.4 

-1933.9 

-142.7 

-171.5 

-121.0 

1.0 
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Table  E-II  Relative  Vorticity  and  Circulation  Components  at  Radial  II 


r' 

w0 

SSI 

^n 

1 

CS2 

drl 

/dr 

dr2/dr 

dr3/dr 

(sec"1) 

(sec"1) 

(sec""1) 

Isec-1) 

(ft/sec) 

(ft/sec) 

(ft/sec) 

.3 

-520 

+697.3 

-89.6 

+756.6 

+55.8 

-4 

.26 

-3.0 

-520 

.35 

-520 
-520 

.4 

-520 

I 

I 

i 

i 

i 

n 

-624.1 

+782.3 

-29.9 

+802.3 

+59.2 

-1.4 

-1.0 

.45 

-780.1 

+910.3 

+59.3 

+870.8 

+64.2 

+2.8 

+1.9 

-656.9 

+809.4 

-11.1 

+816.7 

+60.2 

-.5 

-.3 

.5 

-729.2 

+869.2 

+30.6 

+848.8 

+62.6 

+1.4 

+1.0 

-802.9 

+929.0 

+72.4 

+880.8 

+65.0 

+3.4 

+2.4 

.55 

-767.3 

+899.8 

+52.0 

+865.2 

+63.8 

+2.4 

+1.7 

-826.3 

+948.2 

+85.8 

+891.1 

+65.7 

+4.0 

+2.8 

.6 

-758.9 

+893.0 

+47.2 

+861.5 

+63.5 

+2.2 

+1.6 

-860.1 

+975.9 

+105.2 

+905.9 

+66.8 

+5.0 

+3.5 

.65 

-895.7 

+1005.0 

+125.6 

+921.5 

+68.0 

+5.9 

+4.2 

-1030.1 

+1115.1 

+202.6 

+980.6 

+72.3 

+9.6 

+6.8 

.7 

-1011.6 

+1100.0 

+192.0 

+972.4 

+71.7 

+9.1 

+6.4 

-1716.1 

+1677.0 

+595.7 

+1281.7 

+94.5 

+28.3 

+20.0 

.75 

-2840.7 

+2598.2 

+1240.1 

+1775.4 

+131.0 

+59.0 

+41.7 

-3509.2 

+3145.7 

+1623.1 

+2068.9 

+152.6 

+77.2 

+54.6 

.8 

-3738.6 

+3334.5 

+1755.1 

+2170.0 

+160.1 

+83.5 

+59.1 

-2283.6 

+2141.8 

+920.9 

+1530.9 

+112.9 

+43.8 

+31.0 

.85 

-741.6 

+878.8 

+37.3 

+853.9 

+63.0 

+1.7 

+1.2 

0 

+271.4 

-387.6 

+528.4 

+38.9 

-18.4 

-13.0 

.9 

+741.6 

-336.0 

-812.5 

+202.8 

+14.9 

-38.6 

-27.3 

+1522.4 

-975.5 

-1259.9 

-139.9 

-10.3 

-59.9 

-42.4 

.95 

+5608.9 

-4322.8 

-3601.4 

-1933.9 

-142.7 

-171.4 

-121.3 

+7179.8 

-5609.5 

-4501.6 

-2623.5 

-193.6 

-214.2 

-151.7 

1.0 
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Table  E- III  Relative  Vorticity  and  Circulation  Components  at  Radial  III 


r' 

uf 

i  . 

csi. 

^ 

il, 

?S2, 

dr-j/dr 

dr2/dr 

dr3 

/dr 

(sec"' ) 

(sec"1 ) 

(sec" ') 

(sec"1) 

(ft/sec) 

(ft/sec) 

(ft/sec) 

.3 

+2080.3 

-1432.5 

-1579.6 

-384.8 

-28.4 

-75.1 

-53.2 

.35 
.4 

.45 

.5 

1 

1 

i 

j 

' 

■ 

' 

i 

1 

0 

+271.4 

-387.6 

+528.4 

+38.9 

-18.4 

-13.0 

.55 

-3120.5 

+2827.4 

+1400.4 

+1898.2 

+140.0 

+66.6 

+47.1 

-4160.7 

+3679.4 

+1996.4 

+2354.9 

+173.7 

+95.0 

+67.2 

.6 

-5200.8 

+4531.3 

+2592.4 

+2811.5 

+207.4 

+123.3 

+87.3 

.65 

1 

' 

1 

! 

i 

' 

1 

i 

' 

\ 

1 

i 

.7 

-3649.9 

+3261.0 

+1703.7 

+2130.7 

+157.2 

+81.0 

+57.4 

-4478.9 

+3940.0 

+2178.8 

+2494.6 

+184.1 

+103.7 

+73.4 

.75 

-5149.8 

+4489.6 

+2563.2 

+2789.1 

+205.8 

+122.0 

+86.3 

-5681.4 

+4925.0 

+2867.8 

+3022.5 

+223.0 

+136.5 

+96.6 

.8 

-5380.8 

+4678.8 

+2695.5 

+2800.5 

+213.3 

+128.3 

+90.8 

-2588.0 

+2391.2 

+1095.3 

+1664.5 

+122.8 

+52.1 

+36.9 

.85 

0 

+271.4 

-387.6 

+528.4 

+38.9 

-18.4 

-13.0 

+423.2 

-75.2 

-630.0 

+342.6 

+25.2 

-29.9 

-21.3 

.9 

+1112.4 

-639.7 

-1025.0 

+40.0 

+2.9 

-48.7 

-34.5 

+2283.6 

-1599.0 

-1696.1 

-474.1 

-34.9 

-80.7 

-57.1 

.95 

+7018.4 

-5477.3 

-4408.1 

-2552.6 

-188.3 

-209.8 

-148.5 

+9090.3 

-7174.4 

-5596.3 

-3462.2 

-255.5 

-266.3 

-188.5 

1.0 

- 
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